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ABSTRACT 


The  development  of  the  turbulent  compressible  boundary  layer  on  two 
typical  helicopter  rotors,  for  a  range  of  hover  conditions,  has  been  cal¬ 
culated  using  two  different  analytical  methods.  The  techniques  used  were 
the  Differential  Method,  which  employs  the  differential  form  of  the 
boundary  layer  momentum  equations  and  solves  for  the  local  velocity 
gradients  before  integrating  them  over  the  surface  to  follow  the  flow 
development,  and  the  Integral  Method,  which  employs  the  integrated  form 
of  the  momentum  equations,  having  assumed  a  velocity  profile  form,  and 
solves  for  the  development  of  the  characteristic  boundary  layer  thickness 
parameters  and  skew  angle.  Both  methods  decouple  the  chordwise  and  span- 
wise  boundary  layer  equations  without  making  any  small  crossflow  assump¬ 
tions.  The  effects  of  rotational  speed,  vortex  induced  crossflows, 
surface  curvature  and  applied  chordwise  pressure  gradients  were  evaluated 
separately  and  in  combination  to  simulate  rotor  airfoil  boundary  layer 
growth.  The  effect  of  the  rotation  was  found  to  be  small.  In  all  but 
the  pressure  gradient  cases,  boundary  layer  development  along  the  stream- 
wise  direction  followed  closely  two-dimensional  behavior.  The  results 
from  the  Differential  Method  indicated  that  the  influence  of  unfavorable 
chordwise  pressure  gradients  is  reduced  in  rotating  flows  by  the  presence 
of  an  inward  spanwise  velocity  component  aft  of  the  quarter  chord  line. 
The  influence  of  typical  airfoil  upper  surface  curvature  was  found  to  be 
negligible.  Correlation  between  the  results  predicted  by  the  two  methods 
is  good  for  the  characteristic  thickness  parameters  and  the  surface  shear 
stresses  but  poor  for  the  calculated  skew  angles.  The  use  of  either,  or 
both,  of  these  methods  should  improve  the  prediction  of  surface  skin 
friction  effects  and  give  a  more  reliable  definition  of  separation  point 
in  a  rotating  flow. 
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1.  INTRODUCTION 


With  the  notable  exception  of  von  Karmen' s  exact  solution  for  the  flow 
over  a  rotating  disk.  Reference  (l),  all  the  early  attempts  at  solutions 
of  the  boundary  layer  form  of  the  Navier-Stokes  equations  were  restricted 
to  flows  that  could  be  generally  classified  as  two  dimensional;  that  is, 
flows  where  the  mean  velocity  within  the  boundary  layer  is  in  the  same 
direction  as  the  external  flow.  This  neglect  of  flows  which  make  up  the 
majority  of  those  occurring  in  practice  is  not  surprising,  since  the  number 
of  velocity  gradient  and  shear  stress  terms  which  must  be  considered  in 
the  full  three-dimensional  solution  is  large  and  at  first  sight  intimi¬ 
dating.  Fortunately,  two-dimensional  solutions  were  found  to  be  adequate 
for  most  engineering  applications,  and  activity  was  concentrated  in  this 
area,  with  several  exact  and  empirical  solutions  being  found.  Solutions 
were  based  for  the  most  part  on  the  results  of  the  experimental  programs 
which  paralleled  the  analytical  studies. 

With  the  arrival  of  higher  aircraft  speeds  and  the  adoption  of  3wept  lift¬ 
ing  surfaces,  it  became  evident  that  two-dimens ional  boundary  layer  solu¬ 
tions  would  no  longer  be  ad»quate.  This  led  to  a  widening  of  the  iield  of 
stucty  to  include  the  whole  range  of  realistic  cases  from  the  early  work  on 
the  boundary  layers  on  yawed  cylinders  of  Reference  (2)  to  the  up-to-date 
techniques  discussed  in  References  (3),  (M»  and  (5)  which  make  full  use  of 
modern  computing  machinery  and  allow  the  calculation  of  boundary  layer 
development  not  only  over  simple  surfaces  but  also  around  involved  shapes. 
Paralleling  the  development  of  techniques  to  calculate  boundary  layer 
growth  over  yawed  surfaces  was  the  study  of  the  boundary  ^ayer  on  com¬ 
pressor  and  turbine  blades  and  on  aircraft  propellers  am.  helicopter 
rotors.  In  both  the  swept  and  rotating  surface  problem?,  the  three- 
dimensionality  is  introduced  by  the  presence  of  a  spam.ise,  or  radial, 
pressure  distribution  and,  in  the  rotating  case,  by  the  extra  presence  of 
centrifugal  and  Coriolis  accelerations. 

The  early  rotating  disk  work  by  von  K arm an ,  Reference  (l),  and  the  later 
effort  in  the  same  area  in  Reference  (6),  did  little  to  ease  the  task  of 
trying  to  establish  the  behavior  of  viscou3  flow  on  a  rotating  propeller 
or  rotor  blade.  The  disk  flow  is  dominated  by  a  considerable  radial  com¬ 
ponent  and  there  is  no  circumferential  growth  of  the  viscous  region, 
whereas  on  the  blade,  with  its  radially  disposed  leading  edge,  the  chord- 
wise  or  tangential  effects  predominate  and  radial  flow  effects  are  com¬ 
paratively  small. 

The  earliest  significant  study  of  a  rotating  blade  as  opposed  to  a  disk  is 
given  in  References  (7)  and  (8),  and  the  analogous  problem  of  flow  on  the 
floor  of  a  curved  channel  is  given  in  Reference  (9)«  Dater  work  in 
Reference  (10)  investigated  the  special  case  of  flew  over  a  helical  blade, 
and  most  recently  work  in  References  (ll),  (12),  and  (13)  has  discussed 
the  development  of  boundary  layers  on  helicopter  rotors  in  forward  flight. 

The  basic  problem  that  all  of  these  investigators  have  had  to  face,  in 
addition  to  the  nonlinearity  of  the  governing  equations ,  is  that  the  pair 
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of  boundary  layer  momentum  equations  used  to  define  the  flow  are  coupled, 
and  the  assumptions  made  to  decouple  the  equations  that  the  methods 
differ.  The  approach  used  in  Reference  (8)  was  to  consider  velocities 
and  gradients  of  velocity  in  the  radial  direction  to  be  of  small  order 
compared  to  chordwise  or  tangential  velocities  and  gradients.  This  allows 
the  chordwise  equation  to  be  decoupled  from  the  spanwise  equation,  and, 
for  laminar  flows  at  least,  a  modified  Blasius  solution  may  be  developed 
in  the  chordwise  direction.  The  limitation  of  this  approach  is  that  it 
may  only  be  applied  on  blades  of  high  aspect  ratio  sued  at  large  distances 
from  the  center  of  rotation.  Later  workers,  particularly  Reference  (9). 
assumed  some  characteristic  velocity  profile  form  and,  on  providing  some 
coupling  between  the  streamwise  and  normal  velocity  components,  solved 
the  equations  of  motion  using  integral  techniques.  In  References  (ll) 
and  (13),  the  small  crossflow  work  of  Reference  (9)  is  used  as  the  de¬ 
parture  point  in  a  regular  perturbation  expansion  which  allows  the  calcu¬ 
lation  of  the  development  of  the  boundary  layer  velocity  profiles  in  cases 
where  the  translational  effects  are  superimposed  on  the  steady  rotation 
and  enables  the  effects  of  unfavorable  pressure  gradients  to  be  included. 

A  Blasius  solution  for  the  primary  or  chordwise  flow  is  still  used, 
however. 

The  historical  background  discussed  above  is  not  intended  to  be  all- 
inclusive,  but  merely  to  be  indicative  of  the  principal  approaches  to  the 
problem  of  calculating  the  development  of  the  boundary  layer  on  a  rotating 
blade.  These  solution  techniques,  with  the  exception  of  Reference  (9), 
are  restricted  to  laminar  flows.  It  is  felt  by  the  authors  that  this  may 
not  be  a  realistic  state  of  affairs  since  for  most  full-scale  helicopter 
applications,  the  Reynolds  number  of  the  flow  is  well  above  that  required 
for  transition,  even  on  a  smooth  flat  plate,  and  it  is  considerably 
greater  than  that  for  transition  on  the  aerodynami cally  rough  airfoils  in 
practical  use.  Another  restriction  is  that  all  the  presently  available 
techniques  deal  only  with  incompressible  flows,  an  inappropriate  assumption 
in  light  of  the  high  subsonic  Mach  numbers  found  on  present-day  helicopter 
rotor  blades.  Consequently,  it  was  felt  that  if  new  methods  of  calculation 
were  to  be  set  up  or  if  existing  methods  were  to  be  developed  further,  the 
effects  of  the  turbulent  nature  of  the  flow  and  compressibility  should  be 
taken  into  account.  In  addition,  it  was  felt  that  the  unique  capabilities 
of  modern  computing  machines  should  be  exploited  as  far  as  possible. 

Two-dimensional  boundary  layer  theory  has  advanced  step  by  step  with  ex¬ 
perimental  investigation  of  actual  flows.  Unfortunately,  no  significant 
body  of  experimental  data  is  yet  available  to  the  worker  in  three-dimen¬ 
sional  flows,  and  consequently,  none  of  the  theories  developed  or  under 
development  can  be  considered  valid  until  experimental  proof  is  found. 

With  this  in  mind  and  the  above  requirements  as  a  gride,  it  was  decided  to 
develop  two  separate  and  independent  methods  which  could  be  used  for  com¬ 
parison  and  corroboration  of  the  results . 

Since  the  method  of  Reference  (9)  was  the  only  one  of  the  early  three- 
dimensional  solutions  to  consider  the  calculation  of  turbulent  flows,  it 
was  decided  to  use  this  as  the  basis  for  the  integral  method.  Modifica¬ 
tions  included  allowing  for  the  effects  of  compressibility  and  adverse 


pressure  gradients.  The  resulting  method  vas  developed  for  application 
with  modem  computing  methods  which  were  unavailable  at  the  time  of  the 
original  work. 

For  the  second  method,  it  was  decided  to  use  an  extension  of  a  technique 
which  has  received  considerable  use  in  two-dimensional  applications.  This 
is  to  carry  out  a  full  solution  of  the  differential  equations  for  the 
chordwise  gradients  of  velocity  and  to  integrate  the  gradients  progress¬ 
ively  over  the  surface.  In  the  three-dimensional  case,  if  the  velocity 
components  are  known  along  a  radial  line,  the  spanwlse  velocity  gradients 
can  be  calculated  and  the  chordwise  and  spanwlse  flow  equations  can  be 
evaluated  simultaneously.  This  differential  method  has  the  potential  to 
provide  a  more  exact  solution  of  the  viscous  flow,  but  it  still  relies  on 
certain  assumptions  concerning  the  shear  stress  distribution  for  which  no 
three-dimensional  turbulence  data  is  available. 

Using  both  the  integral  and  differential  methods,  the  development  of  the 
boundary  layer  over  two  representative  rotors  vas  calculated  for  a  range 
of  conditions  including  various  angular  velocities,  vortex  induced  cross¬ 
flows,  adverse  pressure  gradients,  and  blade  surface  curvatures.  The 
development  and  use  of  these  methods  and  the  results  are  discussed  in  the 
following  sections. 
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2.  DEVELOPMENT  OF  THEORY 


2.1  THE  GOVERNING  EQUATIONS 

2.1.1  The  Momentum  Equation 

Reference  ( il+ )  gives  the  generalized  vector  form  of  the  equation  of  motion 
for  an  incompressible  viscous  fluid  moving  relative  to  a  fixed  reference, 
frame  as 


£ 

—  (pq)  +  JjpVq2  -  pqxw  *  -  Vp  -  yVx (l) 
where  q  is  the  general  velocity  vector  and  w,  the  vorticity,  is  given  by 

w  =  Vxq 

If  the  axis  system  is  now  rotated  with  angular  velocity  fl,  all  ;erms ,  with 
the  exception  of  the  unsteady  contribution,  apply  unchanged  relative  to 
the  rotating  frame.  The  relation  between  rates  of  change  of  a  vector 
transformed  from  a  fixed  to  a  rotating  frame  is  given  in  Reference  (15)  as 

ft  ^pi^i  =  ft  (pSJr  +  2£xPi  +  p2.x(RxE'  +  |t^P5.)x£. 


The  vector  r  is  the  position  vector  of  a  point  in  the  rotating  frame.  If 
only  steady  flow  in  the  rotating  system  is  considered,  then  the  above 
equation  reduces  to 


3t  ^P1^I  =  ^LxpS.  + 


The  first  term  on  the  right-hand  side  represents  the  Coriolis  and  the 
second  the  centrifugal  force.  When  this  is  substituted  into  Equation  (l), 
the  general  equation  of  motion  for  a  fluid  in  a  rotating  axis  system 
becomes 


2flxpq  +  pflx(nxr)  +  JgpVq2  -  pqx^  =  -  Vp  -  pVxu 


(2) 


I  II  III  IV  V  VI 


Again  following  Reference  (ll*),  the  component  equations  of  motion  will  be 
written  relative  to  an  orthogonal  curvilinear  axis  system  with  axes  x,  y 
and  z.  The  scaling  factors  will  be  hi,  h£  and  hj,  which  will  generally  be 
functions  of  x,  y  and  z,  and  the  unit  vectors  i_,  and  k  respectively. 

The  vector  quantities  may  be  written  in  terms  of  their  components, 


r  =  xi_  +  yj_  +  zk , 


k 


_ _  ,  _ 


a  ■  ui  +  Vj_  +  Wk 


2  '  n*l  +  V  +  “ti 


The  vorticlty  may  also  be  expanded  In  like  form  with 

“  =  “xi  +  0)yJ.  +  u>zk 

"°*  VOrtiClt3r  ls  deflMd  “  ^  to  VX,  and  may  also  be  written 

hli  h2j.  h^k 

4  =  — i _  1_  9_  3_ 

hlh2h3  3x  3y  g  z 

h_^u  h2v  h^w 

with  components  uv  =  — i —  [i./.  v  3  .  ,  1 

*  *2>>3  [  37  ^3W  *  a?  0>2»)  j 

[  fc  (hlu)  *  fc  <»3»)  ] 

“2  "  hlh2  ffe  <h2T)  '  J7  <>'1u)  j 

^r^^e^'on^Lr  by  tOT  lDt°  ltS  P-U. 

2t(wQy  _  VflZ)i_ 

+  2p(uflz  -  vflx)j_  ( 

+  2P(vflx  -  ufiy)k  , 

Term  II,  pfix^xr),  gives 

p  +  zat)ax  -  (0^3  ♦  nz2)x]i  I 

+  o  [<*»x  ♦  «»,)«,  -  (V  ♦  n^y  JA  j  (l 

♦  0  [<*>*  *  ya,)az  -  (0,2  ♦  o^ja  j  £ 

Since  the  operator  7  in  curvilinear  coordinates  is  8iven  by 

7  =  i-i_i  +  l_i_,  l_3 

hl  3x  -  h2  3y  1  +  h3  37  £ 
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and  q2  =  u2  +  v2  +  v2 ,  term  III  may  be  expanded  to  give 

■  hi  [“il  *  V  *  w£]  i  I 

o  I  3vi  3v  A  3v"l  .  ( 

+  r—  ur—  +  vr—  *■  v—  J  \ 
h2  l  3y  3y  3yJ  *  / 

■  p  r  3u  3v  .  3vl  .  I 

h3  L  3z  3z  3zJ  —  I 


.  3v  . 

+  vr—  +  v 
3z 


Taking  the  component  form  of  the  vorticity  u  and  expanding  it,  together 
with  the  components  of  the  velocity  q,  gives  term  IV, 

Pi*—  *  P  (VU)2  -  vti>y )i_  ) 

*  p(vux  -  uwz)j_  /  (6 

♦  P(lhi)y  -  Vlux)k  ’ 

The  pressure  term,  term  V,  is  simply 

Vp  =  +  L_lEj[  (7 

p  hi  3x  -  h2  3y  1  h3  3z  -  u 


with  the  viscous  effects,  as  represented  by  term  VI,  being  given  by 


[*  (h3“z)  -  ll  (h2«-y^ 

i  \ 

[k (hiu,x)  ■  k  (h3wz) 

L  > 

[k  •  k (hlUx). 

*  / 

hlh2 


Collecting  the  and  k  groups  from  the  expanded  terms  gives  the  com¬ 
ponent  equations  of  motion  in  the  x,  y  and  z  directions  respectively; 


x  Equation: 


2p[vfly  -  vflz]  ♦  p  [(yOy  ♦  zflz)nx  -  (fly2  +  flz2)x]  -r  [u|^  +  vg+  vf“] 

-  p  [w*  *  wuy]  "  -  h7^‘  hfe[l7<h3^)-  fc<h2^)]  C9x) 
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y  Equation: 

2p  |unz  -  vflxj  +  p  pxflx  +  zflz)ny  -  (ftz2  + 

-  “  [w“x  -  H  '  '  9^  ‘ 


z  Equation: 

2p  |vflx  -  uflyj  +  p  ^(xftx  +  yfly)nz  -  (nx2  + 
-  „  [u«v  -  v%]  "  -  i  If  ♦ 


At  this  point  in  the  development  of  two-dimensional  boundary  layer 
equations,  it  vould  .ormal  to  carry  out  an  order  of  magnitude  compar¬ 
ison  of  terms  and  rev.  .ce  the  equations  to  the  standard  short  form. 

Before  this  is  done  in  the  present  case,  it  is  useful  to  consider  the 
surface  over  which  the  flow  is  taking  place,  the  orientation  of  the  axes 
relative  to  the  surface,  and  the  dependence  of  the  scaling  factors  on  x, 
y  and  z.  Figure  (l)  gives  the  system  of  axes  used  in  the  present  work. 
Following  References  (8)  and  (ll),  and  others,  the  axes  will  be  fixed  in 
the  rotating  surface.  The  x  and  y  coordinates  will  define  the  surface 
z  =  0,  and  z  will  follow  the  local  normal  to  the  surface.  Curvature 
will  be  restricted  to  the  x  axis,  making  scaling  factors  h_  and  h_equal 
unity.  Since  the  blade  may  be  represented  by  a  cylindrljal  surface, 
the  scale  factor  h^  becomes  a  function  of  z  alone. 

It  is  conventional,  in  treating  the  boundary  layer  equations  in  two 
dimensions  with  the  flow  in  the  x  direction,  to  assume  that  x,  u  and  p 
sure  all  of  order  unity,  0  (l),  sued  that  z  and  w  axe  of  the  order  of  the 
boundary  layer  thickness,  0(6).  In  addition  to  these,  the  existence  of 
the  crossflow  field  introduces  terms  involving  y,  v,  hi,  h2»  h3  and  Q(, 
which  are  all  of  order  0(l),  and  fly,  and  ftx  or  order  0(6).  The  angular 
velocity  ftz  is  generally  large  and  will  be  treated  as  equal  to  ft,  the 
angular  velocity  of  the  rotor  about  its  axis.  By  expanding  equations  9x, 
9y  and  9z  and  discarding  terms  on  the  basis  of  their  relative  orders  of 
magnitude,  and  at  the  same  time  substituting  ftfor  ftz  and  replacing  the 
scaling  factors  h2  and  h3  with  their  value  of  1.0,  the  full  equations  of 
motion  sure  reduced  to  the  equations  of  motion  for  the  steady  flow  of  a 
viscous  fluid  over  a  singly  curved  rotating  surx'ace: 


7 


x  Equation: 
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This  completes  the  development  of  the  boundary  layer  momentum  equations. 
2.1.2  The  Continuity  Equation 

The  generalized  vector  form  of  the  equation  of  continuity  in  a  fixed  axis 
frame  is  given  by 


+  7*(pq)  =  0  (11) 

Since  the  density  p  is  not  a  vector  quantity,  Equation  (ll)  will  apply 
unchanged  if  the  coordinate  system  is  rotated.  If  the  further  requirement 
is  made  that  the  flow  relative  to  the  rotating  axis  system  not  be  time- 
dependent,  then  the  continuity  equation  for  steady  flow  relative  to  a 
rotating  axis  system  is  simply 


V-(pq)  =  0  (12) 

Bearing  in  mind  that  in  the  curvilinear  coordinate  system  the  operator  7 
is  given  by 
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_1  _3 
hi  3x 
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Equation  (12)  may  be  expanded  into 


V*(p<l)  =  (l^hgpu)  +  ~  (h3hipv)  +  ~  (hih2Pw)  =  0  (13) 

In  defining  the  surface  and  coordinate  system  for  the  rotor  blade  problem, 
hj  and  h3  were  specified  above  as  unity  and  hi  as  a  function  of  z  alone. 
With  this  in  mind  the  continuity  equation  (Equation  ( 13 ) )  now  becomes 
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2.1.3  The  Energy  Equation 


No  original  treatment  of  the  energy  equation  will  be  attempted  here;  the 
solution  of  the  equation  derived  in  Reference  (l4)  will  be  used.  For  the 
case  of  the  adiabatic  wall  to  be  assumed  here,  a  relation  for  the  local 
temperature  is  given  of  the  form 


T  =  T  + 


i  r  - 

_±_  u* 
2Cp  |_  ° 


V2  - 


(15) 


If  this  is  rearranged  to  give  the  local  density  in  terms  of  the  local 
velocity  and  the  Mach  number  of  the  external  flow,  it  becomes 


(16) 


2.1.4  Reynolds  Stress  Effects 

The  flow  equations  developed  above  have  been  in  terms  of  the  instantaneous 
values  of  the  flow  variables.  To  determine  the  effect  of  the  turbulent 
velocity  fluctuations  in  the  boundary  layer,  it  is  conventional  to  replace 
the  instantaneous  values  of  the  flow  variables  with  the  sum  of  the  time 
mean  and  fluctuating  parts.  For  example, 

u  a  u  ♦  u' 


The  bar  denotes  a  time-mean,  and  the  dash  a  fluctuating  quantity.  The 
mechanics  of  performing  this  operation  are  covered  in  References  (l6)  and 
(17),  and  will  not  be  dwelt  on  in  this  analysis.  In  the  absence  of  any 
significant  body  of  experimental  turbulence  data  in  three-dimensional 
flows,  only  those  terms  known  to  be  significant  in  two-dimensional  flows 
will  be  considered,  together  with  their  obvious  cross-flow  partners. 
Bearing  this  in  mind  and  neglecting  triple  and  higher  order  products  of 
the  fluctuating  terms ,  the  motion  and  continuity  equations  become 
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The  most  significant  effect  of  turbulence  is  the  addition  of  an  apparent 
shear  stress  term  tp  each  of  the  momentum  equations.  These  are  the 
Reynolds  stress  terms,  the  time  mean  of  the  product  of  the  fluctuating 
velocities  in  the  x  and  y  directions,  and  the  normal  velocicy  fluctuation. 
The  precise  form  of  these  terms  and  their  dependence  on  the  mean  flow 
parameters  vill  be  discussed  in  a  later  section. 

2.2  THE  SOLUTIONS 

2.2.1  Limiting  Conditions 

The  equations  governing  the  flow  over  a  cylindrical  surface  rotating  about 
some  axis  perpendicular  to  its  length  were  developed  in  Equations  (l7x) 
and  (l7y)t  the  momentun  equations,  in  Equation  (l8),  the  continuity 
equation,  and  in  Equation  (16),  a  solution  of  the  energy  equation  for  the 
special  case  of  adiabatic  flov.  In  succeeding  sections,  they  vill  be 
manipulated  and  subsequently  solved  to  determine  the  development  of  the 
boundary  layer  over  the  surface.  The  two  methods  of  solution  used  differ 
considerably,  but  they  must  both  satisfy  the  same  boundary  conditions: 
that  there  be  zero  slip  at  the  vail,  z  ■  0,  vith  u  ■  v  ■  0  and  consequently 


3u  3u  3v  3v  _ 

——  b  —  b  ■—  a  —  a  0 

3x  3y  3x  3y 

In  the  case  of  the  flov  over  a  solid  vail  vith  no  transpiration,  v0  is  also 
equal  to  zero.  A  further  condition  is  that  at  large  distances  from  the 
surface,  u  and  v  must  approach  U*  and  V*,  the  velocities  of  the  external 
flov. 


2.2.2  The  Differential  Method 


It  is  convenient  at  the  outset  to  nondimens ionalize  the  velocity  components 
in  the  x  (chordvise)  and  y  (spanvise)  directions  by  dividing  by  U,,,  the 
velocity  of  the  external  flov  in  the  chordvise  direction.  They  become 
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The  equations  of  motion  written  in  terms  of  c  and  s  become 
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(In  the  interest  of  clarity,  the  bars  have  been  dropped  from  the  time  mean 
quantities  except  for  the  (pw  +  p  'v' )  term. )  The  viscous  and  Reynolds 
stress  terms  have  been  combined  in  the  terms  3tx/3z  and  3iy/3z  above. 

The  main  object  of  the  solution  is  to  isolate  the  gradients  of  velocity  in 
the  chordvise  direction,  3c/3x  and  3s/3x,  and  these  are  the  principal  un¬ 
knowns  in  the  above  equations.  Before  they  can  be  determined,  however, 
substitutions  must  be  made  for  the  pressure  gradient  terms,  the  normal 
velocity  term  (fly  +  pTwr),  and  the  shear  stress  variation  normal  to  the 
surface. 

The  pressure  gradient  terms  may  be  found  by  evaluating  Equations  (l7x)  and 
(I7y)  at  the  surface  in  the  absence  of  viscous  effects  or  by  considering 
the  differentiated  form  of  Bernoulli's  equation  in  a  rotating  system. 
Bernoulli '8  equation  for  the  flew  being  considered  here  may  be  written  as 
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If  this  is  differentiated  with  respect  to  x  and  y,  it  is  found  that 
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(20) 


The  term  containing  the  normal  velocity  may  be  removed  by  substitution 
from  Equation  (l8),  the  continuity  equation.  Equation  (l8)  must  first  be 
rearranged  and  integrated  along  a  normal  to  the  surface  to  give 
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“  (pu)  +  hi  (fv)  j  dz  (21) 
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Equation  (16)  gives  an  expression  for  the  density  in  terms  of  the  velocity 
components,  and  if  this  is  differentiated  with  respect  to  x  and  y,  the  inte¬ 
grand  of  Equation  (21)  becomes 
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Collecting  terms  involving  the  spanvise  and  chordwise  velocity  gradients 
and  expressing  the  velocities  in  terms  of  the  velocity  ratios  c  and  s 
leads  to  the  final  form  of  Equation  (21): 

(57  ♦  rzn-  -  ^  fj  [<£)  ♦  c2Fj]  ♦  £  C.F2  +  f*  MhjF2 
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The  shear  stress  terms  in  Equations  (l9x)  and  (l9y)  must  now  he  expressed 
in  terms  of  the  mean  flow  quantities.  From  Equations  (l7x)  and  (lTy),  the 
shear  stress  terms  3tx/3z  and  3iy/3z  are  given  by 


and 


Knowing  the  velocity  profiles,  the  viscous  shear  stresses  may  be  determined 
directly.  Unfortunately,  some  assumptions  must  be  made  to  obtain  an 
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expression  for  the  Reynolds  stress  tens.  In  two-dimensional  boundary 
layer  solutions,  the  relation  between  turbulent  shear  stress  and  mean  flow 
is  most  commonly  supplied  by  the  eddy  viscosity  or  related  mixing  length 
hypotheses.  Other  approaches  involve  solutions  of  the  turbulent  kinetic 
energy  equation  or  an  eddy  viscosity  transport  equation.  However,  as  a 
result  of  the  present  poor  state  of  knowledge  of  the  behavior  of  turbu¬ 
lence,  all  methods  employ  empirical  relationships  to  some  degree. 

Reference  (18)  gives  an  up-to-date  summary  of  all  the  methods  in  common 
use  and  an  assessment  of  their  applicability.  In  three-dimensional  flows, 
the  data  available  is  very  limited  and  the  relationship  between  the  tur¬ 
bulent  shear  stress  and  the  mean  flow  is  still  a  matter  of  argument . 
Considering  this,  it  will  be  assumed  for  the  present  solution  that  a  mix¬ 
ing  length  relationship  may  be  applied  and  the  turbulent  shear  stress  may 
be  related  to  the  normal  gradient  of  the  local  re  a  ul  tan’  /cloclty  by 

where  r  *  (u2  +  v2)1  Is  the  local  resultant  velocity  and  pr'wk  is  the 
resultant  turbulent  shear  stress.  The  two  do  not  necessarily  act  in  the 
same  direction.  Reference  (  5  )  suggests  that  it  is  most  probable  that  the 
turbulent  shear  stress  acts  in  the  direction  of  the  maximta  rate  of  strain 
of  the  mean  flow  given  by 


(25) 


With  the  magnitude  of  the  reaultant  shear  stress  given  by  Equation  (2l») 
and  the  direction  by  Equation  (25),  it  may  be  further  asaumed  that 
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It  is  now  possible  to  arrive  at  expressions  for  the  Reynolds  stress  com- 


ponents  in  the  x  and  y  directions: 
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The  form  of  the  mixing  length  £  suggested  for  two-dimensional  flow  in 
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R^ferenTe  (19)  bu  be«n  used  In  Um  present  work.  It  is  given  by 


y  •  lull  iff.) 


where  im  la  th«  vtior  »r  1  In  the  outer  boundary  layer,  typically  lm/6  * 
0.09  to  0.095,  and  *  la  the  von  Karman  conatant  relating  the  mixing  length 
to  dlatance  normal  to  the  wall.  The  two-dimensional  value  of  «  is 
typically  0.k,and  this  will  be  used  here.  The  relations  given  in 
Equations  (26)  only  apply  In  that  part  of  the  boundary  layer  where  tur¬ 
bulent  stresses  predominate.  In  the  regions  close  to  the  wall,  viscous 
effects  control  the  flow,  and  the  expression  for  the  nixing  length  has  to 
be  modified.  This  has  been  done  for  two-dimensional  flows  in  Reference 
(20),  but  It  was  felt  that,  because  of  the  approxlmat  ‘one  made  in  the  tur¬ 
bulent  region  and  the  three-dimensionality  of  the  flow,  the  refinement 
suggested  would  be  Inappropriate.  Instead,  the  shear  stress  grsdients  in 
the  laminar  region  were  determined  by  Interpolating  between  conditions  at 
the  wall  and  at  the  edge  of  the  sublayer. 


Since  flow  In  this  region  Is  controlled  by  the  effects  of  viscosity,  and 
References  (20  sod  (27)  bl£illght  the  si^ilflcance  of  the  shear  velocity, 
u*.  close  to  the  wall,  this  was  used  as  the  basis  for  the  Interpolation. 
The  shear  stress  gradients  are  given  by 
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The  Interpolation  is  given  In  full  In  the  appendix. 

The  expressions  for  the  pressure  gradient.  Equation  (20),  the  term 
( PV  ♦  o *w ' ) ,  Equation  (.'?),  and  the  shear  stress  gradient.  Equation  (26), 
may  now  be  substituted  Into  the  boundary  layer  momentum  equations,  (I9x) 
and  ( 19y ) ,  to  give  the  unknewn  chordwlse  gradients  in  terms  of  known 
parameters: 
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The  shear  stress  gradients  3xx/3z  and  3Ty/3z  have  been 
form  since  Equations  (26)  are  strictly  only  applicable 
laminar  sublayer.  Equations  (27x)  and  (27y)  represent 
with  two  unknowns,  and  as  such  they  are  soluble.  They 
through  the  integral  term  which  contains  both  unknowns 
equations.  The  group  p^U^2  may  be  removed  as  a  factor 
rearranged  to  give 
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Equation  (31)  may  be  expressed  in  terms  of  a  new  unknown, 

G(z)=bit£+b2f 

°(z)  =  (bjf^  +  b2fyl)  J G(z)dz  +  (bxfx2  +  b2fy2)  (32) 
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This  may  be  easily  solved  for  0  using  iterative  techniques  which  are 
ideally  suited  for  machine  computation.  The  problem  has  now  been  reduced 
to  the  solution  of  two  simultaneous  equations  in  3 c/3x  and  3s/ 3x: 

fyl  If  -  fxl  ft  =  fx2fyl  “  fy2fxl  (33) 

bl|f  +  b2  |f  =  G(z)  (3k) 


The  solution  is  now  a  routine  matter  of  separating  the  unknowns  using 
Equations  (33)  and  (31*). 

The  calculation  of  the  gradients  3c/3x  and  3s/3x  described  above  is 
carried  out  at  a  number  of  points  along  the  normal  to  the  surface  within 
the  boundary  layer.  The  actual  number  of  solution  points  is  determined  by 
the  need  to  accurately  specify  the  velocity  profiles,  and  their  spacing  is 
determined  by  the  need  to  adequately  describe  the  rates  of  change  of 
quantities  normal  to  the  wall.  In  the  present  calculation,  the  points  have 
been  distributed  logarithmically  from  the  wall  with  a  minimum  of  twenty 
points  within  the  boundary  layer.  It  is  then  possible  to  integrate  the 
chordwise  gradients  to  the  next  point  using  a  forward  stepping  finite 
difference  technique.  If  this  id  done  for  a  number  of  stations  along  a 
constant  x  line,  for  which  the  flow  conditions  are  known,  then  the  develop¬ 
ment  of  the  boundary  layer  over  the  surface  may  be  calculated. 

2.2.3  The  Integral  Method 

The  principal  difficulty  in  turbulence  work  lieB  in  specifying  the  shear 
stress.  Most  work  thus  relies  on  some  type  of  semiempirical  model  to 
define  the  turbulent  shear  stress  throughout  the  boundary  layer.  By  using 
an  integral  technique,  this  necessary  approximation  may  be  partially 
avoided.  The  momentum  equations  may  be  integrated  with  respect  to  the 
distance  normal  to  the  surface  out  to  the  free  stream.  The  integral  of  the 
shear  stress  term  is  the  value  of  the  skin  friction  at  the  surface  for 
which  an  empirical  model  must  then  be  provided.  This  can  be  done  without 
affecting  the  accuracy  of  the  results  in  many  cases,  and  at  the  same  time, 
numerical  complexity  can  be  significantly  reduced.  However,  variations 
normal  to  the  surface  are  eliminated  because  the  boundary  layer  is  defined 
in  teims  of  overall  parameters.  This  technique  has  proven  its  usefulness 
in  two-dimensional  work,  and  its  extension  to  three  dimensions  is  shown  to 
be  practical  in  Reference  (9). 

Several  steps  must  be  taken  before  integration  of  the  momentum  equations 
becomes  feasible.  The  pressure  gradient  terms  can  be  evaluated  directly 
from  Equations  (lTx)  and  (l7y)  by  inserting  the  free  stream  density  and 
velocities.  The  effect  of  variation  of  pressure  normal  to  the  surface  is 
considered  to  be  negligible.  Thus,  the  resulting  pressure  gradient  ex¬ 
pressions  are  applicable  from  the  free  stream  to  the  surface.  The  value  of 
the  shear  stress  term  outside  the  boundary  layer  is  zero.  The  pressure 
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gradient  expressions 
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When  these  expressions  are  substituted  into  the  momentum  equations  and 
written  with  the  intended  integration,  the  following  equations  result: 
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Since  the  external  flow  is  considered  to  be  parallel  to  the  surface,  terms 
involving  Wm  have  been  set  equal  to  zero.  The  bars  have  been 
dropped  from  the  time  mean  quantities  except  for  the  (pw  +  p'w')  term.  In 
the  x  component,  the  two  terms  with  coefficients  of  (pw  +  p  'w' )  have  been 
combined. 


Equation  (lO),_the  continuity  equation,  may  be  employed  to  supply  an  ex¬ 
pression  for  (pw  +  p'w')  in  terms  of  the  remaining  variables,  as  was  done 
in  section  2.2.1.  The  result,  Equation  (21),  is  repeated  here. 

z 

(pw  +  p'w' )  =  -  J  £  (pu)  +  hi  |^-  (pv)j  dz  (21) 

*0 

Now  (pw  +  p'w')  may  be  eliminated  from  the  momentum  equations,  and  the 
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resulting  double  integral  terms 
6 

-  [  (&±pD  Mpi)  r  a^] 

0  hl  92  •*  hl  3z  J  L  3*  hl  3y  F  dz  (37x) 

0  0  J 

and  „ 

r  _  - 6  -  \ 

-  (pw  f  |1  d2  =  /  3v  r  a^u)  8(0U)1 

J  3z  J  L3x  +hig^Jdzdz  (3 

may  be  integrated  by  parts  in  the  same  manner  as  their  ^  < 

counterparts.  The  actual  change  in  the  sc^le  fLtnl  l  two‘di“«f ional 
equivalent  to  a  typical  boundary  layer  is  In  i  l  °Ver  a  di8tance 
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conUenu?ty8equatdion"0f  sfnce^he^sheL^t8'0^0  integral  tennS  COMB  from  the 

stream,  their  integrals  have in  the  free 
the  chordwise  and  spanwise  skin  frirtioL  t  the  lower  llmits‘  These  are 
can  be  rewritten  as  ’  T°x  and  Toy’  These  equations 
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n!?8*  Uof8  convenient  to  transfonn  the  distance  noimal  to  the 
all  by  using  a  Stevartson  transformation  of  the  form 
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The  integral  terms  can  then  be  written  in  an  incompressible  form 
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where  A  is  the  thickness  of  seme  related  incompressible  boundary  layer. 
These  are  the  three-dimensional  forms  of  the  displacement  and  momentum 
thicknesses.  The  terms  6*  and  0*x  are  'the  same  as  in  two-dimensional 
boundary  layer  theory.  The  term  6p  represents  a  change  in  thickness  due 
to  compressibility.  The  other  thickness  parameters,  6y,  0*y»  8yy  and  0jX» 
arise  from  the  existence  of  the  velocities  in  the  y  direction  and  are  not 
present  in  two-dimensional  flows.  If  the  density  ratio  can  be  related  to 
the  local  velocities,  the  only  unknowns  remaining  on  the  left  side  of  the 
momentum  equations  are  the  integral  terms  involving  u  and  v.  The  solution 
of  the  energy  equation  (Equation  ( 16) )  can  be  used  to  supply  this  ratio. 


In  two-dimensional  work,  only  the  momentum  and  displacement  thicknesses  in 
the  x  direction  appear.  A  shape  factor,  H  =  5*/0xx»  can  be  used  to  relate 
these  variables,  thereby  permitting  a  solution  of  the  momentum  equations 
without  a  knowledge  of  the  velocity  profiles.  In  the  three-dimensional 
case,  velocity  profiles  are  required  to  relate  the  numerous  thickness 
parameters  and  reduce  the  number  of  variables  in  the  integrated  equations. 
Power  law  profiles  of  a  form  suggested  in  Reference  (  q  ), 


and 


where 


-A*  (A) 
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1/N 


02} 


2  1/N 

r=  (1-  A}  (N-  *  * 
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=  — -  tan( e )  * 


•ox 
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axe  used  in  the  streamwise  and  crossflow  directions.  Their  shapes  are 
shown  for  various  exponents  in  Figure  (  2  )  •  The  angle  e  is  the  angle 
between  the  surface  streamline  and  the  external  flow.  Since  the  overall 
thickness  parameters  rather  than  the  details  of  the  velocity  profiles  are 
required  in  the  equations,  these  profiles  are  considered  to  be  sufficiently 
representative  without  involving  the  complexity  of  integrating  a  three- 
dimensioned  form  of  a  perhaps  more  rigorous  _uw  of  the  wedl/law  of  the 
wedee  profile.  Thus,  power  law  profiles  axe  resolved  into  streeunwise  and 
chordwise  directions. 
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and  are  then  used  to  evaluate  the  integrals  in  terms  of  the  thickness  A, 
the  skew  angle  e,  and  the  quantity  N: 
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The  momentum  equations,  normalized  by  p^lt^2,  now  appear  as 
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Two  things  remain  to  be  specified:  the  skin  friction  and  the  quantity  N 
from  the  power  law.  In  zero  pressure  gradient  cases,  a  value  of  seven  is 
commonly  used  for  N;  however,  for  more  general  use,  1/N  may  be  related  to 
the  shape  factor  H  by  evaluating  6*/Qxx  in  two-dimensional  form,  giving 


1  _  H-l 
N  "  2 


(vr) 


The  change  in  H  as  the  boundary  layer  develops  in  pressure  gradient  cases 
can  be  obtained  from  an  empirical  relation  such  as  that  given  by 
von  Doenhoff  and  Tetervin,  Reference  ( 21 ) » 
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The  momentum  thickness,  6q,  is  evaluated  using  the  streamwise  velocity 
profile.  There  are  several  other  methods  of  accounting  for  shape  factor 
changes,  but  the  above  relationship  is  the  most  readily  compatible  with 
the  present  integral  technique.  Since  the  conditions  under  which  it  was 
determined  are  similar  to  the  cases  with  which  it  will  be  used,  good 
results  might  be  expected.  These  conditions  include  turbulent  boundary 
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layers  under  moderate  adverse  pressure  gradients  on  the  upper  surface  of 
an  airfoil. 


The  skin  friction  must  also  he  obtained  from  an  empirical  relation. 
Integral  techniques  have  most  commonly  relied  on  the  Ludweig-Tillman  law 
from  Reference  (22)» 


.123  x  10 


-0.678H 
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This  relationship  illustrates  the  dependence  of  the  skin  friction  on  the 
shape  factor,  H.  However,  since  the  work  of  von  Doenhoff  and  Tetervin 
was  based  on  the  skin  friction  law  of  Squire  and  Young,  Reference  (23), 
their  relation 
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must  be  used.  Since  H  does  not  enter  into  this  relation,  it  may  be  in¬ 
accurate  for  large  values  of  the  shape  factor.  However,  the  growth  of  the 
calculated  boundary  layer  should  be  unaffected  by  this  since  it  is  con¬ 
trolled  by  the  empirical  shape  factor  relationship. 

Both  the  shape  factor  relation  and  the  skin  friction  law  are  two- 
dimensional  in  origin  and  are  to  be  applied  along  the  streamwise  direction. 
The  skin  friction  can  be  resolved  into  spanwise  and  chordwise  components. 


This  follows  from  the  definition  of  $.  The  shape  factor  is  not  a  direc¬ 
tional  property.  Its  value  is  considered  to  be  unique  at  every  point  and 
therefore  determinable  from  the  two-dimensional  relation. 

A  solution  of  the  momentum  equations  is  now  possible.  This  becomes  clear 
if  they  are  rewritten  as  two  simultaneous  partial  differential  equations 
in  the  form 
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If  Initial  values  of  A,  4,  and  H  are  specified  along  a  lino  parallel  to 
the  leading  edge,  the  spanwlse  gradients  can  be  evaluated  and  the  skin 
friction  and  shape  factor  gradient  can  be  calculated  at  regular  points 
along  this  line.  The  moaentia  equations  can  then  be  solved  for  the  chord- 
vise  gradients  at  each  of  these  points;  these  gradients  nay  then  be  inte¬ 
grated  in  the  x  direction  from  each  point  to  determine  the  parameters  A, 

♦  *f>d  H  along  a  new  line  parallel  to  the  leading  edge.  The  process  can 
then  be  repeated  as  far  downstream  as  required.  The  value  of  6p  can  be 
used  to  transform  the  results  back  into  a  compressible  coordinate  system. 

2.3  THE  BOUHDAHY  COW  DITTOES 


2.3.1  The  External  Flow 


The  two  Independent  methods  of  predicting  turbulent  boundary  layer  develop¬ 
ment  on  a  rotor  blade  .equire  that  the  external  flow  magnitude  and  direc¬ 
tion  be  specified  beforehand.  This  flow  is  comprised  of  the  rotational 
velocity,  the  velocity  Induced  by  the  pressure  distribution  on  the  blade, 
and  a  superimposed  crossflow  velocity  component  due  to  the  close  passage 
of  the  trailing  tip  vortex  from  the  preceding  blade. 

In  a  rotating  rectangular  coordinate  system,  the  spanvise  and  chordvise 
rotational  velocity  coq>onents  can  be  expressed  as 

Ua  «=  fly  (51a) 

and 

Va  ^  (51b) 

respectively;  y  is  measured  from  the  axis  of  rotation,  as  described  in 
Figure  (1).  In  performing  the  calculation  on  a  rotor  blade,  the  x  *  c/1* 
line  is  made  to  coincide  with  the  quarter  chord  line.  Spanvise  rotational 
flow  it  outward  ahead  of  this  line  (negative  x)  and  inward  behind  this 
line  (positive  x). 

The  flow  over  a  rotating  flat  plate  is  described  completely  by  the  above 
rotational  velocities.  When  the  surface  assumes  a  thickness  distribution 
such  as  that  of  an  airfoil,  the  local  flow  velocities  are  altered  by  the 
resulting  pressure  distribution.  The  pressure  distribution  is  expressed 
in  terms  of  the  pressure  coefficient, 


Cp  = 


which,  for  later  use,  can  be  rewritten  as 


1  + 


*5  P„Ua2cP 
Pa 


C52) 


26 


The  subscript  (a)  refers  to  the  ambient  conditions,  which  are  also  the 
static  conditions  in  the  rotating  system,  and  the  subscript  (<■)  refers  to 
the  local  conditions  just  outside  the  boundary  layer.  In  the  absence  of 
exact  pressure  distribution  data  which  includes  the  effect  of  the  external 
crossflows,  an  approximate  two-dimensional  pressure  distribution  based  on 
the  chordwise  flow  is  used. 

To  solve  for  the  local  velocity,  U^,  in  terms  of  the  pressure  coefficient 
and  the  ambient  conditions,  the  pressures  sure  first  related  to  the  Mach 
numbers , 
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This  may  be  combined  with  Equation  (52)  to  eliminate  the  pressure  p-t 
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The  local  Mach  number  may  then  be  expressed  as 
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The  local  velocity  can  be  obtained  from  *  M^a^,  but  first  the  local 
sound  speed,  a  ,  must  be  defined.  The  adiabatic  temperature  ratio  and 
Equation  (53)  are  used  to  find  an  expression  for  the  ratio  of  local  to 
ambient  sound  speeds, 
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The  local  velocity  is  then  written  in  terms  of  the  pressure  coefficient 
and  the  static  conditions, 
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By  using  the  adiabatic  density  ratio,  the  local  density,  pm, 
obtained  in  a  similar  manner  and  expressed  as 
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1  + 


**  paUa2Cp 


can  be 


(56) 


The  local,  velocity  and  density  are  thus  known  in  terms  of  the  ambient  con¬ 
ditions  and  the  pressure  distribution  of  the  airfoil.  Typical  upper  sur¬ 
face  pressure  distributions  at  angles  of  attack  of  0°  and  U°  for  a  Mach 
number  of  0.65  from  Reference  (2I*)  are  plotted  in  Figure  (3).  They  were 
obtained  from  two-dimensional  tests  conducted  on  a  12jf-thick  symmetrical 
airfoil  section.  In  both  cases,  an  approximation  of  the  measured  pressure 
distribution  has  been  made  by  assuming  a  constant  pressure  over  the  first 
quarter  of  the  chord  and  an  increasing  pressure  over  the  remainder.  Since 
the  actual  pressure  distributions  are  also  functions  of  the  Mach  number, 
the  approximations  ere  meant  to  be  representative  of  the  adverse  pressure 
gradients  that  may  bt>  encountered  on  an  airfoil.  These  pressure  gradients 
cause  a  deceleration  of  the  flow  and  a  thickening  of  the  boundary  layer 
toward  the  trailing  edge.  It  is  this  region  of  an  airfoil  cn  which  a 
turbulent  boundary  layer  is  most  likely  to  separate.  Lower  surface 
pressure  distributions  at  positive  angles  of  attack  may  be  adequately  pre¬ 
sented  by  a  constant  or  mildly  adverse  pressure  distribution. 

The  magnitude  of  the  pressure  gradient  may  be  obtained  from  the  slope  of 
the  pressure  coefficient,  and  if  this  is  nondimeusionalized  by  the  chord, 
the  approximated  distributions  have  values  of  dCp/d(x/C)  of  1  and  2.  The 
local  velocity  gradient  to  which  the  boundary  layer  is  subjected  is  ob¬ 
tained  by  differentiating  the  velocity  expression.  Equation  (55),  to  give 
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Similarly,  the  density  gradient  is  given  by 
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In  multibladtd  rotors,  a  significant  external  spanwise  flow  may  be  imposed 
by  the  presence  of  a  trailing  tip  vortex  from  the  preceding  blade.  It  has 
been  shown  in  References  (25)  and  (26)  that  such  vortices  remain  near  the 
plane  of  the  rotor  until  they  peas  the  following  blade  and  are  then  swept 
down  away  from  the  blades.  The  region  of  blade-vortex  interference  is 
somewhat  inboard  of  the  tip  due  to  wake  contraction;  for  this  analysis,  it 
can  be  assumed  to  be  centered  at  the  90-percent  radial  station.  A  more 
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precise  positioning  of  the  vortex  would  be  necessary  for  rotor  performance 
calculations . 

The  magnitude  of  this  spanwise  flow  can  be  estimated  from  potential  theory 
by  writing  the  stream  function  of  a  line  vortex  and  its  image,  each  of 
which  is  parallel  to  the  x  direction  at  a  distance  Z  from  the  surface 
z  =  0  and  at  a  distance  Y  from  y  =  0, 


T  (Y-y)2  +  (Z-z)2 

*  =  ^  Se  (y-y)2  +  (Z+z)2 


(59) 


The  resultant  spanwise  velocity  at  the  surface  is  determined  by  differen¬ 
tiating  the  stream  function  with  respect  to  z  and  setting  z  equal  to  zero. 


v  _  li  _  1 1  z2 

00  ~  3z  "  TT  Z  (Y-y)2  +  z2 


(60) 


The  portion  of  the  vortex  which  is  in  the  region  of  the  blade  is  repre¬ 
sented  as  a  straight  segment  parallel  to  the  x  direction.  This  spanwise 
velocity  distribution  may  be  superimposed  on  the  chordwise  flow  to  create 
a  representative  rotor  blade  flow  to  which  a  boundary  layer  may  be  sub¬ 
jected.  The  thickness  of  the  boundary  layer  on  the  surface  has  been 
neglected  in  f  .  i  acing  V*,  since  it  is  much  less  than  the  height  Z,  which 
itself  has  on.  ■v  ■  .  approximate  value. 

2.3.2  The  Rotor  Blade 

In  the  development  of  the  equations  of  motion  and  of  the  two  solution 
methods,  the  scale  factor  depicting  coordinate  curvature  in  the  chordwise 
direction,  h^,  was  retained.  This  scale  factor  is  related  to  the  radius 
of  curvature  of  the  x  coordinate  in  the  z  direction  by 


h1  3z 


Since  h,  has  been  taken  to  be  solely  a  function  of  z,  the  derivative  can 
be  written  as  a  toted  rather  them  a  partial,  and  the  scede  factor  can  be 
determined  as  a  function  of  z  by  integration.  This  gives 
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Rn  ♦  Z 


(61) 


where  R„  has  been  written  as  the  radius  of  curvature  of  the  surface  plus 
the  distance  from  the  surface,  R0  +  z,  and  the  constant  of  integration 
has  been  determined  by  letting  ^  *  1  at  i  =  0.  In  the  integral  method, 
one-half  of  the  boundary  layer  thickness,  6,  is  used  as  the  appropriate 
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value  of  z.  The  gradient  of  the  scale  factor  can  be  obtained  by 
differentiation. 


dhi  _i_ 

dz  *  R0 


(62) 


and  can  be  seen  to  be  constant  above  every  point  on  the  surface.  If  the 
surface  is  assumed  to  be  a  circular  arc,  R0  will  be  fixed  for  all  points 
on  the  blade.  Such  an  arc  is  shown  in  Figure  (4),  along  with  an  NACA  0012 
airfoil  for  comparison.  The  two  curves  are  tangent  at  the  30-percent 
station,  and  they  intersect  at  the  trailing  edge.  It  can  be  shown  by  trig¬ 
onometry  that  the  radius  of  such  an  arc  is  R0  =  4. 1C,  where  C  is  the  chord 
in  feet.  This  approximation  will  be  used  in  conjunction  with  the  assumed 
pressure  distribution  and  the  vortex  crossflow  representation  to  simulate 
actual  conditions  for  the  prediction  of  turbulent  boundary  layer  growth  on 
a  rotating  blade. 
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3.  APPLICATION  OF  THEORY 


3.1  DEVELOPMENT  OF  THE  NUMERICAL  METHODS 


3.1.1  The  Differential  Method 

With  the  exception  of  some  difficulty  with  the  sheen:  stress  gradients  the 
development  of  the  numerical  methods  for  the  differential  technique  was 
comparatively  routine.  A  computer  program  was  written  first  to  solve  the 
reduced  form  of  the  equations  of  motion,  Equations  (33)  and  (3^+),  for  the 
local  chordwise  velocity  gradients  along  normals  to  the  surface  in  the 
initial  X  plane,  and  then  to  integrate  them  over  the  surface  to  follow  the 
flow  development. 

A  simple  fixed  mesh  structure  was  used  as  a  basis  for  the  solution.  The 
spacing  of  the  solution  points  in  the  direction  of  solution,  the  X  or 
chordwise  direction,  wets  determined  by  a  compromise  between  several 
factors:  the  sensitivity  of  the  integration  procedure  to  step 
size;  the  requirement,  noted  in  Reference  (5  )»  that  for  a  stable  solution 
the  ratio  of  the  chordwise  aid  spanwise  step3,  Ax/Ay,  should  not  exceed 
Ua/Voo;  and  the  more  practical  economic  factor  that  excessively  small  chord- 
wise  steps  would  result  in  unacceptably  long  computation  times.  In  the 
early  stages  of  development,  a  fixed  chordwise  step  size  of  0.5  times  the 
smallest  boundary  layer  thickness  present  on  the  surface  was  used.  It  was 
found  that  this  could  be  increased  considerably,  and  an  initial  step  size 
of  0.01  ft  was  eventually  used  for  all  the  conditions  examined.  This 
initial  step  size  was  used  for  the  first  four  steps  of  the  calculation  to 
allow  conditions  to  stabilize  and  was  then  increased  by  a  factor  of  2*  for 
the  remainder  of  the  calculation.  The  regular  step  size  was  felt  to  be 
acceptable  in  view  of  the  simple  forward  stepping  integration  technique 
being  used.  If  an  iterative  technique  were  employed,  or  if  an  implicit 
solution  were  being  considered,  then  variable  chordwise  steps  would  pex'haps 
be  more  appropriate.  The  step  size  in  the  spanwise  direction  was  also 
chosen  to  give  economic  computation  times  without  sacrificing  definition  of 
the  flow  field  in  regions  of  high  spanwise  rates  of  change  or  compromising 
the  stability  requirement  noted  above.  The  steps  used  were  1.0  ft  on  the 
large  blade,  Ay/R  =  0.0?5,  and  0.5  ft  on  the  small  blade.  Ay /R  =  0.05.  In 
a  differential  solution,  and  particularly  in  one  in  which  numerical  dif¬ 
ferentiation  of  the  variation  of  parameters  normal  to  the  surface  has  to 
be  carried  out,  the  distribution  of  points  normal  to  the  surface  must  be 
considered  with  particular  care.  Foints  must  be  distributed  with  not  only 
the  rates  of  change  of  conditions  normal  to  the  wall  in  mind  but  also  the 
peculiarities  of  the  differentiating  techniques  to  be  used.  In  order  to 
keep  the  computation  as  simple  as  possible,  a  fixed  normal  coordinate  system 
was  used  with  the  points  distributed  exponent! ally  away  from  the  wall 
following 

^  .  0.5xl0-'*[l00-1(k-1)-1] 

The  dimensions  of  z  axe  in  ft,  and  k  is  an  integer  with  value  1  at  the  wall 
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and  UO  at  the  point  furthest  from  the  vail.  This  distribution  gave  at 
least  20  points  within  the  boundary  layer  for  the  thinnest  case  considered, 
with  a  minimum  of  8  points  in  the  laminar  sublayer. 

Evaluation  of  the  coefficients  of  Equations  (33)  and  (3^)  requires  that 
the  gradients  of  several  quantities,  in  a  normal  and  spanwise  direction, 
be  known.  Since  the  mesh  size  is  regular  in  the  spanwise  direction  and 
rates  of  change  are  comparatively  low,  a  simple  three-point  Lagrangian 
interpolation  formula  was  used  to  determine  the  gradients  in  this  direc¬ 
tion.  The  nonuniform  step  size  normal  to  the  wall  presented  some 
problems,  however.  Initially,  a  five-point  Lagrangian  interpolation 
scheme  with  nonuniform  intervals  was  used.  It  was  found  that  this  method 
introduced  significant  errors  in  the  first  and  unacceptable  errors  in  the 
second  derivative  because  of  the  concentration  of  points  in  the  lower  end 
of  the  groups.  The  problem  was  overcome  and  the  errors  were  minimized  by 
transforming  the  normal  coordinate  into  a  form  where  the  intervals  between 
solution  points  were  uniform  and  carring  out  the  differentiation  using  a 
second  order,  orthogonal  Lagrangian  polynomial  fit  on  five  adjacent  points. 
The  transformation  used  was  simply  the  inverse  of  the  expression  used 
above  to  define  zk,  the  differentiation  being  carried  out  in  terms  of  k. 

In  order  to  give  a  basis  for  comparison,  the  two  calcualtion  methods  used 
identical  initial  conditions.  The  boundary  layer  was  fully  developed  and 
turbulent  from  the  leading  edge,  and  initially,  both  the  skin  friction 
and  the  total  boundeuy  layer  thickness  were  given  by  the  simple  zero 
pressure  gradient  forms 
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Since  the  external  flow  is  specified,  the  above  expressions  may  be  used  to 
give  the  skin  friction  and  boundary  layer  thickness  along  the  starting 
line,  x/c  ■  0.1.  The  law  of  the  wall/law  of  the  wake  form  for  the  turbu¬ 
lent  profile  suggested  in  Reference  (27)  and  the  tanh  form  of  the  sublayer 
profile  suggested  in  Reference  (28)  were  used,  together  with  the  surface 
skin  friction  and  boundary  layer  thickness  given  above,  to  define  the 
initial  velocity  profile.  The  initial  crossflow  velocity  profile  was 
assumed  to  be  given  by  Equation  (1*3).  The  calculation  was  found  to  be 
insensitive  to  the  initial  value  of  the  surface  skew,  ^  in  Equation  (1*3), 
the  spanwise  velocity  profile  being  dominated  by  the  spanwise  ccmponent  of 
the  8treamvl8e  profile,  and  an  initial  value  of  zero  was  subsequently  used. 

The  most  significant  effort  during  the  development  of  the  method  was  in 
finding  satisfactory  form  for  the  shear  stress  gradient  normal  to  the 
wall.  It  was  originally  intended  that  the  sheeir  stress  gradients  in  the 
turbulent  part  of  the  layer  be  determined  by  differentiating  the  stresses 
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calculated  using  the  mixing  length  approach,  Equation  (26).  These  would 
then  be  faired  into  the  wall  values  using  an  interpolation  based  on  the 
tanh  series  suggested  in  Reference  (19)  and  extended  to  the  three-dimen¬ 
sional  case.  This  method  was  found  to  be  impractical  since  it  was  very 
sensitive  to  the  quality  of  the  match  at  the  edge  of  „ne  sublayer.  Any 
nonuni formity  in  the  shear  stress  gradient  normal  to  the  feeds  into  the 
calculation  of  the  chordwise  velocity  gradient,  and  this  tends  to  magnify 
the  irregularity  as  the  solution  proceeds.  This  is  associated  with  the 
difficulty  of  obtaining  precise  gradients  using  a  numerical  differentiation 
scheme  discussed  above.  The  problem  is  accentuated  because  obtaining  the 
gradient  of  the  shear  stress  involves  differentiating  the  product  of 
several  previous  differentiation  operations,  Equations  (26).  The  errors 
are  compounded  and  the  difficulty  of  obtaining  a  uniform  match  at  the  edge 
of  the  sublayer  arises.  Considerable  effort  was  put  into  finding  a 
solution  to  this  problem,  but  with  little  success.  Eventually,  the  first 
five  points  outside  the  sublayer  were  fitted  by  a  least  mean  square 
process  to  follow  u+2,  and  the  interpolation  outlined  in  the  appendix  was 
used. 

3.1.2  The  Integral  Method 

The  development  of  the  integral  method  entailed  no  major  problems.  A 
computer  program  was  initially  written  to  evaluate  the  integral  form  of 
the  momentum  equations  (Equations  (50x)  and  (50y))  simultaneously  over  a 
simulated  rotor  blade  surface.  The  Ludweig-Tillman  skin  friction  relation 
was  originally  used  in  this  program  to  calculate  boundary  layer  growth 
with  zero  pressure  gradient.  The  shape  factor  was  assumed  to  have  a  con¬ 
stant  value  of  H  *  1.286.  This  corresponds  to  a  power  law  profile  with  an 
exponent  of  1/7,  as  indicated  by  Equation  (1*7).  The  required  initial 
values  of  6  were  obtained  from  a  two-dimensional  relationship  suggested 
in  Reference  (l6). 

6  .  0.3TX  [Sa]  5 

This  presupp  see  a  turbulent  boundary  layer  from  the  leading  edge  (x  *  0) 
onward,  as  was  previously  indicated. 

The  initial  value  ox  the  tangent  of  the  skew  angle,  $,  is  not  as  readily 
estimated.  At  first,  this  was  done  by  assuming  $  *  0.0  along  the  first 
row  of  chordwise  stations.  This  appeared  reasonable ,  since  the  solutions 
that  followed  over  the  remainder  of  the  blade  indicated  that  the  crossflow 
stayed  small.  Nevertheless,  several  initial  values  for  $,  including  nega¬ 
tive  numbers ,  were  tried;  the  solutions  were  found  to  >e  unaffected  after 
four  or  five  steps.  Thus,  ^  =  0.0  has  been  used  in  starting  all  sub¬ 
sequent  calculations. 

With  these  initial  conditions  and  a  knowledge  of  the  external  flow,  the 
integral  momentum  equations  "an  be  evaluated  according  to  the  procedure 
outlined  at  the  end  of  Section  2.2.3.  No  spanwise  boundary  conditions  are 
required;  all  variables  are  deterrrir  ~d  either  from  the  preceding  chordwise 
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station  or  from  the  external  flow  conditions.  The  spanwise  gradients  are 
calculated  with  a  three-point  Lagrangian  interpolation.  The  end  gradients 
are  determined  from  the  last  three  spanwise  points  with  an  appropriate 
form  of  the  interpolation,  which  is  Justifiable  except  when  the  spanwise 
flov  is  entering  the  calculation  region,  as  pointed  out  in  Reference  (5). 
This  condition  occurs  at  the  tip  aft  of  the  quarter  chord  due  to  the 
circular  form  of  the  streamlines.  The  possible  ill  effects  of  this  are 
considered  to  be  small  because  the  external  flow  and  its  gradients  are  well 
defined  in  these  regions  and  because  they  behave  smoothly  in  all  cases. 

The  next  major  step  in  the  development  of  the  integral  method  was  the  use 
of  a  variable  shape  factor  to  allow  consideration  of  adverse  pressure 
gradients.  Equations  (55)  through  (58)  are  used  to  specify  the  free  stream 
velocity  and  density  and  their  chordwise  gradients  when  pressure  distri¬ 
butions  of  the  form  discussed  in  Section  2.3.1  are  applied.  The  computa¬ 
tion  procedure  was  modified  to  calculate  the  shape  factor  gradient  from  the 
von  Doenhoff-Tetervin  relation  before  solving  the  motion  equations.  The 
shape  factor  at  each  step  is  determined  from  its  value  and  gradient  at  the 
previous  step.  In  constant  pressure  regions,  the  gradient  is  zero.  The 
shin  friction  law  was  replaced  with  the  Squire-Young  relation,  Equation 
(49),  as  is  required  for  compatibility  with  the  von  Doenhoff-Tetervin 
equation.  The  change  of  skin  friction  relationships  decreased  4  in  zero 
pressure  gradient  cases.  This  is  examined  more  fully  in  the  discussion  of 
results. 

The  chordwise  step  size  was  chosen  to  be  0.1  foot.  Steps  which  were  an 
order  of  magnitude  smaller  were  tried,  but  no  significant  difference  was 
observed  without  pressure  gradients.  For  adverse  pressure  gradient  cases, 
the  O.i-foot  step  is  larger  than  von  Doenhoff  and  Tetervin  recommend,  but 
the  l  ocation  of  the  separation  point  on  a  blade  with  a  chord  of  2  feet  is 
changed  less  than  5  percent.  Larger  chordwise  steps  would  significantly 
impair  the  results.  Spanwise  gradients  were  found  to  be  very  small  in 
relation  to  chordwise  gradients.  Thus,  the  use  of  a  larger  spanwise  step 
size  is  considered  appropriate.  A  step  of  1  foot  was  used  in  this 
direction . 

3.2.1  Dat'i  Input 

Data  input,  to  the  boundary  layer  programs  is  accomplished  via  punched 
cards.  This  allows  the  greatest  ease  and  flexibility  in  setting  various 
conditions  for  the  computation.  Details  of  the  usage  of  the  programs  are 
covered  in  the  Documentation  of  Software  Report.  The  data  required  by  Doth 
programs  is  basically  the  same.  Ambient  air  conditions,  including  density, 
temperature,  pressure,  and  sound  speed,  are  read  in.  The  blade  geometry  is 
specified  by  inputting  the  blade  radius,  chord  and  the  blade  surface  radius 
of  curvature,  along  with  the  inboard  radial  station  for  starting  the  cal¬ 
culation.  Spanwise  and  chordwise  step  sizes  are  also  read  in.  The  angular 
velocity  combined  with  the  minimum  pressure  coefficient  and  the  nondimen- 
sionalized  pressure  gradient  enable  the  programs  to  calculate  the  chordwise 
external  flow  field.  A  nonrotating  case  can  be  run  by  setting  fl  »  0.0  and 
giving  a  positive  value  for  the  two-dimensional  free-stream  velocity. 

(Use  of  both  this  velocity  and  a  nonzero  fi  will  not  simulate  forward 
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flight.)  The  amount  of  the  blade  which  is  subjected  to  constant  pressure 
in  the  adverse  pressure  gradient  cases  is  an  input  variable,  but  values 
other  than  one  quarter  of  the  chord  have  not  been  used.  The  vortex  cross- 
flow  is  defined  by  supplying  values  of  the  circulation  strength  and  the 
spanwise  and  normal  positions  of  the  vortex  core.  In  addition,  the 
integral  method  requires  the  initial  value  of  the  shape  factor.  The  data 
output  format  is  selected  in  each  method  by  a  series  of  inputs  which  are 
discussed  in  the  Documentation  of  Software  Report. 

A  wide  range  of  input  data  is  acceptable,  although  a  few  reservations 
should  be  observed.  Singularities  that  should  be  avoided  occur  when  the 
external  chordwise  velocity  is  zero  and  when  the  surface  radius  of  curva¬ 
ture  is  zero.  The  first  will  happen  when  both  the  two-dimensional  velocity 
and  0  are  zero.  Free  stream  flows  should  generally  remain  subsonic.  Even 
though  the  solution  of  the  energy  equation  which  has  been  used  is  appro¬ 
priate  above  these  velocities,  no  consideration  has  been  given  to  the 
possibility  of  shock  formation  or  shock  wave/boundary  layer  interaction. 

The  general  choice  of  step  sizes  has  been  discussed  in  Section  3.1. 

In  performing  the  calculations,  standard  sea  level  conditions  were  used  for 
the  ambient  air.  Since  the  effects  of  rotation  are  to  be  considered  rela¬ 
tive  to  helicopter  applications,  two  blade  radii,  UO  feet  and  10  feet,  were 
selected  as  representative  of  large  present-day  helicopter  main  and  tail 
rotors.  The  large  and  small  blades  have  chords  of  2  feet  and  1  foot, 
respectively.  Angular  velocities  of  10,  15,  and  20  rad/sec  for  the  large 
blade  and  Uo,  60,  and  80  rad/sec  for  the  small  blade  were  chosen.  These 
speeds  encompass  a  wide  range  of  angular  velocities  without  introducing 
locally  supersonic  f 1  jw  even  in  the  strong  pressure  gradient  cases.  The 
form  of  the  pressure  distribution,  along  with  the  magnitudes  of  the  surface 
radii  of  curvature,  was  discussed  in  Section  2.3. 

Circulation  strengths  of  + 200  and  +400  ft2/sec  and  a  passage  distance  from 
the  surface  of  Z  ■  2  feet  are  considered  appropriate  for  the  vortex  inter¬ 
fering  with  the  UO-foot  blade.  The  positive  values  represent  a  vortex 
passing  over  the  top  surface  by  imposing  an  outward  flow  and  the  negative 
values  represent  a  vortex  passing  under  the  lower  surface  by  imposing  an 
inward  flov.  The  ci  'culation  strengths  for  the  10-foot  blade  are  +25  and 
+50  ft^/sec  at  a  dist  nee  from  the  surface  of  Z  *  1  foot. 

3.2.2  Data  Output 

Calculated  data  is  printed  out  in  tabular  form  for  each  of  the  programs  ,  as 
is  described  in  detail  in  the  Documentation  of  Software  Report.  Data  is 
listed  according  to  chordwi6e  and  spanwise  position.  The  variables  avail¬ 
able  for  output  include  the  various  thickr*ss  related  parameters ,  the  shear 
stresses,  the  skew  angles,  and  the  calculated  external  flow.  The  differ¬ 
ential  method  inherently  has  ix>re  output  information  available  in  the  form 
of  a  listing  of  the  actual  velocity  profiles,  as  well  as  a  breakdown  norma1, 
to  the  surface  of  the  forces  that  influence  the  flow.  The  available  in¬ 
formation  is  used  as  a  practical  aid  In  the  data  analysis.  To  obtain  a 
record  of  the  output  that  can  be  put  to  future  use  by  the  computer,  the  most 
vital  data  was  punched  on  cards.  This  is  detedled  in  the  Documentation  of 
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Software  Report.  This  form  of  output  was  the  source  for  the  data  summary 
presented  in  Tables  I  through  XII  and  the  plotted  data  in  Figures  (5) 
through  (1*0).  In  both  cases,  the  original  output  cards  were  used  to 
supply  the  required  information  which  was  repunched  in  the  necessary 
format.  The  tables  wore  generated  by  listing  these  punched  output  cards. 
The  output  for  the  figures  was  used  to  generate  curves  for  a  numerically 
controlled  drafting  machine.  The  displacement  thickness  curves  in  Figures 
(5)  through  (l6)  and  the  skew  angle  curves  in  Figures  (17)  through  (28) 
were  obtained  by  fitting  a  curve  with  second  derivative  continuity  through 
the  data  points .  Irregularities  in  the  velocity  profiles  near  the  surface 
where  the  calculation  points  are  extremely  close  together  prevented  this 
curve  fit  routine  from  being  used,  so  Figures  (29)  through  (4),  the 
velocity  profiles,  were  drawn  with  straight  line  segments  connecting  each 
point.  The  figures  are  grouped  to  ill:istrate  separately  the  effects  of 
variation  of  angular  velocity,  vortex  crossflow,  pressure  gradient,  and 
surface  curvature.  The  order  of  these  groups  is  repeated  for  both  the 
large  and  small  blades,  and  this  pattern  is  repeated  for  the  presentation 
of  the  displacement  thickness,  the  skew  angles,  and  the  velocity  profiles 
of  the  flow.  The  summary  data  is  similarly  arranged  in  Tables  I  through 
XXI. 


k.  DISCUSSION  OF  RESULTS 


k.l  EFFECTS  OF  ROTATION 


The  equations  of  motion  for  the  development  of  a  turbulent  boundary  layer 
on  a  rotating  blade  were  obtained  in  Section  2.  The  equation  for  the 
chordwise  flow  is  analogous  to  the  two-dimensional  form  of  the  boundary 
layer  equation  with  extra  forcing  and  inertial  terms  added  by  the  presence 
of  a  crossflow  velocity  component  and  the  fact  that  the  whole  system  is 
rotating  about  axis  with  some  angular  velocity  fl.  The  basic  form  of  the 
boundary  layer  equations  used  in  both  solution  methods  is  given  in 
Equations  (l7x)  and  (l7y)»  and  they  are  repeated  below  for  convenience. 
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5he  terms  -pfi2x  and  -pfl2y  in  Equations  (l7x)  and  (l7y)»  respectively,  are 
components  of  the  acceleration  of  the  rotating  coordinate  system  toward 
the  center  of  rotation.  Except  for  the  origin,  all  points  which  are 
stationary  in  the  rotating  system  have  an  acceleration  in  the  inertial 
frame  of  reference.  In  addition,  a  point  moving  in  the  rotating  system 
has  a  Coriolis  acceleration  perpendi  ular  to  its  direction  of  motion  and 
the  axis  of  rotation.  The  terms  -2pvft  and  +2f>ufl  represent  the  chordwise 
and  spanwiae  rates  of  change  of  momentum  of  a  point  as  it  traverses  a  re¬ 
gion  of  varying  rotational  velocity.  These  two  sets  of  acceleration  terms 
represent  the  direct  effects  of  the  rotation.  Further  indirect  effects 
are  associated  with  the  choice  of  the  coordinate  system.  The  flow  is 
being  defined  relative  to  a  rectangular  coordinate  system  fixed  on  the 
blade  surface  and  rotating  about  an  axis  on  the  blade  quarter  chord  line. 
Because  of  this,  the  spanwise  external  flow  velocity  is  positive  ahead  of 
the  quarter  chord  line  and  negative  aft,  and  an  effective  pressure  gradient 
term  is  introduced  into  the  chordwise  equation  even  when  no  chordwise 
pressure  gredient  is  present.  Equation  (20).  The  spanwise  equation,  (l7y), 
is  of  course  dominated  by  the  radial  gradient  of  the  chordwise  velocity 
associated  with  the  rotation. 

The  effect  of  centrifugal  force  on  the  boundary  layer  is  small  and  only 
becomes  significant  close  to  the  wall  when  compressibility  effects  become 
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large-  This  is  apparent  from  Equations  (27x)  and  (27y)  where  the  centrif¬ 
ugal  contribution  is  given  as  -fi2x( p-p^/Iji) and  -n2y(p-p«)  respectively. 

The  Pa,  terms  are  introduced  from  the  pressure  gradient  substitution.  For 
low  Mach  number  flows,  p-Pa,  and.  the  centrifugal  effect  are  essentially  zero. 
The  influence  of  the  Coriolis  force,  on  the  other  hand,  is  significant, 
particularly  on  the  spanwise  flow.  Since  the  chordwise  velocity  u  iB,  for 
attached  flows  at  least,  positive,  the  effect  of  the  Coriolis  force  will 
be  to  decelerate  the  spanwise  flow.  In  the  chordwise  flow,  the  small  effect 
of  Coriolis  force  may  be  either  accelerative  or  decelerative  depending  on 
whether  the  spanwise  flow  is  outward  or  inward.  Again,  in  the  chordwise 
flow,  the  effect  of  the  "pseudo"  pressure  gradient  term,  Voo3Vao/3x  in 
Equations  (20),  introduced  by  the  external  spanwise  flow, is  generally 
small.  In  rotating  flat  plate  flow,  Va,=-fi(x- c/4)  with  3Va>/3x  *  -  ft,  making 
the  effect  of  the  term  decelerative  ahead  of  the  quarter  chord  line  and 
accelerative  behind.  Since  the  effect  is  proportional  to  ft2,  it  may  be 
expected  to  be  more  significant  on  the  smaller  blades  with  the  higher 
rotational  speeds.  The  equivalent  term  in  the  spanwise  equation, 

UcSUco/ay,  is  dominant  and  controls  the  flow  in  this  direction. 

The  first  noticeable  feature  of  the  boundary  layer  on  the  rotating  surface 
is  the  spanwise  variation  of  its  characteristics  due  to  the  distribution 
of  free  stream  velocity  caused  by  rotation.  This  is  most  readily  observed 
from  the  figures  that  compare  results  for  various  angular  velocities, 
Figures  (5),  (ll),  (17),  (23),  (29),  and  (35),  and  from  Tables  I  and  VII. 
The  increase  in  velocity  with  radius  has  a  thinning  effect  on  the  thickness 
parameters  and  results  in  lower  values  of  the  skin  friction,  just  as  an 
increase  in  free  stream  velocity  affects  two-dimensional  boundary  layer 
growth.  The  boundary  layer  also  becomes  thinner  for  a  given  spanwise 
station  as  the  angular  velocity  is  increased,  as  can  be  seen  by  comparing 
parts  (a),  (b),  and  (c)  of  Figures  (5)  and  (ll). 

When  the  chordwise  momentum  and  displacement  thicknesses  are  compared  with 
two-dimensional  results  obtained  at  the  same  chordwise  velocity,  the 
differences  are  almost  negligible.  There  is  a  very  slight  thinning  in  the 
three-dimensional  case  because  the  resultant  velocity  qJ  *  U»2  ♦  V.2 

along  the  chordwise  line  used  for  comparison  is  large*  -han  that  of  the 
two-dimensional  case.  Further  adjustments  to  U«  in  the  two-dimensional 
case  for  the  sake  of  comparison  would  not  be  meaningful  because  the 
differences  in  results  are  already  beyond  the  limits  of  practical  signif¬ 
icance.  The  behavior  of  the  rotating  boundary  layer  in  terms  of  its  growth 
rate  is  virtually  two  dimensional.  The  effects  of  rotation  of  the  chord- 
wise  grcarth  are  negligible  on  both  blades.  The  terms  which  describe  these 
effects  in  the  x  momentum  equation  are  typically  two  or  three  orders  of 
magnitude  less  than  the  shear  stress  term  which  controls  the  development. 

In  Figures  (Sc)  and  (lie),  the  displacement  thickness  cm  be  seen  to  in¬ 
crease  with  span  rather  than  decrease .  This  occurs  on  the  outer  portion 
of  both  the  large  and  small  blades  for  the  highest  angular  velocity  cases 
and  for  the  pressure  gradient  cases,  but  only  when  the  free  stream  velocity 
exceeds  approximately  600  ft /sec.  Its  origin  may  be  traced  to  the  effect 
of  compressibility  and  the  definition  of  the  displacement  thickness.  From 
the  solution  of  the  energy  equation  (Equation  (l6))»tbe  density  within  the 


boundary  layer  can  be  seen  to  decrease  near  the  surface  as  the  Mach 
number  Increases.  The  displacement  thickness , 
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is  thus  increased  as  the  ratio  p/p„  becomes  less  than  one.  This  thickness 
can  be  defined  as  the  distance  by  vhich  the  surface  would  need  to  be  moved 
if  the  flow  were  invlscid  to  produce  the  same  external  flow  character¬ 
istics.  If  the  density  in  the  boundary  Jayer  is  reduced,  then  the  dis¬ 
placement  thickness  will  correspondingly  increase  if  the  velocity  profile 
is  unchanged. 

The  boundary  layer  thickness,  6,  remains  unaffected  because  the  velocity 
distribution  is  not  significantly  affected.  A  similar  increase  doeu  not 
occur  for  the  momentum  thickness.  The  density  ratio  would  appear  from  the 
definition 
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to  reduce  this  quantity,  but  the  effect  is  not  noticeable.  The  momentum 
thickness  is  the  thickness  of  a  unit  depth  of  the  external  flow  which 
contains  a  momentum  flux  deficiency  of  the  boundary  layer.  It  is  thus 
directly  reiatable  to  the  profile  drag  of  the  surface. 

Values  of  the  thickness  parameters  calculated  using  the  two  different 
methods  are  in  relatively  good  agreement  for  the  flat  plate  case  being 
discussed  here.  The  results  of  the  differential  method  shov  somewhat 
more  sensitivity  to  the  effects  of  Q  than  do  those  from  the  integral 
method. 

The  growth  of  the  skew  single,  e,  is  a  function  of  chordwi^e  and  spanviae 
position.  In  the  integral  method,  the  skew  angle  is  practically  indepen¬ 
dent  of  ft  in  the  zero  pressure  gradient  cases.  In  Figure  (4l),  e  is 
plotted  versus  y  for  a  chordwlse  position  of  1  foot  for  both  the  large  and 
small  blades.  The  intermediate  angular  velocity  for  each  blade  was  used. 
This  curve  shows  the  inverse  relationship  between  c  and  y.  Of  particular 
interest  is  the  close  match  between  the  skew  angle  predicted  by  the 
Integral  method  at  the  tip  of  the  10- foot  blade,  which  is  traveling  at 
600  ft/sec,  and  at  the  y  *  10-ft  position  of  the  40-foot  blade,  which  is 
traveling  at  150  ft/sec.  This  occurs  because  in  the  final  integral  form 
of  the  momentum  equations,  (50x)  and  (50/) •  the  terms  in  the  solution  of 
3e/3x  are  either  constant  in  the  spanvise  direction  or  proportional  to  1/y. 

Skew  angles  predicted  by  the  differential  method  are  generally  larger,  in 
an  absolute  sense,  than  those  predicted  by  the  integral  method.  This  is 
not  surprising,  since  c  is  calculated  by  taking  the  difference  between  the 
external  and  wall  flow  directions  and  the  latter  is  calculated  by  obtaining 
a  ratio  of  the  wall  shear  stresses.  In  the  integral  method  the  shear 
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stress  is  analytically  defined,  whereas  in  the  differential  method  it  is 
calculated  from  the  velocity  profiles.  The  differences  betveen  the  values 
of  e  presented  are  within  those  indicated  by  the  shear  stresses. 

The  chordwlse  skin  friction  coefficient  is  higher  in  the  rotating  cases 
with  the  same  chordwlse  velocity.  Since  the  equations  are  normalized  with 
respect  to  p.U*2,  the  streamwise  skin  friction  coefficient  in  the  integral 
method  is  resolved  into  chordwlse  and  spanwise  components  as 


The  tangent  of  the  skew  angle,  4,  lb  usually  positive  while  the  crossflow, 
V„,  is  negative,  and  thus  the  chordwlse  skin  friction  coefficient  is 
greater  than  TOq/p0.U.2  by  at  least  a  factor  of  00,/Uoo;  tox  1b  not  necessar¬ 
ily  less  than  toq. 

The  spanwise  skin  friction  coefficient  indicates  the  relative  amount  of 
spanwise  flow.  From  the  tabulated  results  it  can  be  seen  that  there  is  a 
difference  of  two  orders  of  magnitude  between  the  spanwise  and  the  chord- 
wise  coefficients.  The  relative  sizes  of  these  two  variables  is  controlled 
by  the  relative  size  of  their  respective  external  flow  velocity  components 
and  the  value  of  the  skew  angle,  as  indicated  in  Equations  (50) .  Both  V® 
and  i  remain  small  in  most  cases .  Values  for  the  flat  plate  skin  friction 
coefficients  calculated  by  the  two  methods  are  in  good  agreement  with  the 
differential  method,  again  indicating  a  greater  sensitivity  to  the  effects 
of  external  crossflow  than  does  the  integral  method.  This  is  reflected  in 
smaller  absolute  values  of  the  spanwise  skin  friction  coefficient  at  the 
inboard  stations,  where  the  ratio  of  spanwise  to  chordwise  external  flow 
velocities  is  greatest. 

Figure  (42)  shows  the  chordwise  skin  friction  coefficient  plotted  against 
Reynolds  number  for  a  wide  range  of  conditions  as  calculated  by  the 
differential  method  on  both  large  and  small  blades  and  the  two-dimensional 
values.  Agreement  is  generally  very  good  except  at  the  high  Reynolds 
number  values  where  compressibility  plays  a  significant  role. 

4.2  Effect  of  Vortex  Crossflow 

The  details  of  the  interference  of  a  tip  vortex  shed  from  a  preceding  blade 
can  be  best  observed  by  studying  the  spanwise  distribution  of  6X*  and  e  at  a 
constant  chordwise  position  over  the  outer  portion  of  the  blade.  This  has 
been  plotted  in  Figure  (43)  for  the  40-foot  blade  for  both  methods.  The 
effects  are  not  shown  as  well  in  the  machine-plotted  figures  due  to  the 
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different  scale.  The  displacement  thickness  increases  on  one  side  and 
decreases  on  the  other  side  of  the  centerline  of  the  vortex,  depending  on 
the  sign  of  the  circulation  strength,  r.  The  same  effect  is  manifest  in 
5  and  8X.  Near  the  centerline  of  the  vortex  these  variables  remain  almost 
unchanged  even  though  the  crojsflov  velocity,  V«o,  is  the  highest  in  this 
position.  However,  the  velocity  gradient  9Va>/3y  has  on  inflection  with  a 
zero  value  at  the  vortex  centerline,  and  this  quantity  influences  the 
chordwise  growth  of  the  boundary  layer  and  causes  the  spanwise  inflection 
in  t.ie  thickness  parameters . 

When  P  is  positive,  the  velocity  V<*>  at  the  centerline  is  outward,  or 
positive,  whereas  the  original  undisturbed  flow  would  have  been  negative. 
Inboard  of  the  centerline  there  is  a  divergence  of  the  external  stream¬ 
lines  where  3V«»/3y  is  positive,  and  the  boundary  layer  is  reduced  in 
thickness.  Outboard,  3V«,/3y  is  negative  and  the  streamlines  tend  to  be 
convergent.  In  this  region,  the  boundary  layer  is  thicker.  The  opposite 
situation  occurs  when  T  is  negative.  The  influence  of  3V«,/3y  is  confirmed 
when  the  results  from  the  T/h  *  +_  200  ft/sec  cases  are  examined.  Here, 
3V«,/3y  is  smaller,  and  as  would  be  expected  the  deviations  from  the  T  =  0 
case  are  reduced  proportionately. 

The  skew  angle  predicted  by  the  integral  method  Increases  inboard  of  the 
vortex  centerline  and  decreases  outboard  of  it,  regardless  of  the  sign  of 
T.  Figure  (43)  shows  that  the  variation  is  larger  for  negative  values  of 
r  than  for  positive.  In  the  integral,  solution  for  e,  terms  with  coeffi¬ 
cients  of  V.«DVoo/3y  have  a  dominant  effect.  The  sign  of  this  expression 
will  always  be  positive  inboard  of  the  vortex  and  negative  outboard,  caus¬ 
ing  the  skew  angle  to  develop  as  it  does. 

The  skew  angle  variation  predicted  by  the  differential  method  is  radically 
different  from  that  discussed  for  the  integral  method.  The  reason  for  this 
difference  is  not  hard  to  find.  The  variation  of  e  in  the  differential 
method,  while  being  influenced  by  V®3Voo/3y  like  the  integral  method,  is 
more  strongly  controlled  by  the  radial  variation  of  the  chordwise  skin 
friction.  Figure  (44)  shows  the  sensitivity  of  cfx  predicted  by  the 
differential  method  at  the  trailing  edge  to  the  presence  of  the  vortex. 

The  integral  method  shows  no  comparable  effect.  The  differential  method 
calculates  skew  angle  by  taking  the  difference  between  external  and  wall 
flow  directions  and  determines  the  latter  from  the  ratio  of  cfX/cfy  at  the 
wall.  This  explains  the  form  of  the  e  variation  when  the  radial  variation 
of  Cfx  and  cjy  under  the  vortex  is  considered. 

4.3  Effect  of  Surface  Curvature 

The  characteristics  of  the  boundary  layer  are  virtually  unaffected  by  the 
surface  curvature  tested  on  both  the  large  and  small  blades.  Careful 
examination  of  the  complete  data  listings  reveals  that  there  is  occasion¬ 
ally  an  insignificant  increase  of  the  order  of  0.1  percent  in  the  thickness 
related  parameters.  There  is  no  detectable  change  in  the  surface  shear 
stress  as  predicted  by  either  method.  The  slight  increase  in  the  growth 
rate  is  a  result  of  the  factor  of  1/hi  in  the  chordwise  velocity  gradient 
term  of  Equation  (I7x),  pu/hi( 3u/3x) .  Since  1/hi  is  less  than  1,  the 


magnitude  of  3u/3x  will  be  increased.  This  term  is  generally  negative; 
the  increased  magnitude  produces  a  thickening  of  the  boundary  layer.  From 
Equation  (6l),  the  value  of  hi  is  seen  to  be  equal  to,  or  less  than, 

(i  +  S/R0).  The  changes  noted  in  6  are  consistent  with  this. 

Similar  effects  were  noted  on  blades  with  vortex  crossflows  and  in  the 
cases  which  combined  surface  curvature  and  adverse  pressure  gradients.  In 
general,  for  the  surface  radii  tested  in  this  report,  the  effects  of  sur¬ 
face  curvature  are  negligible.  The  radius  of  curvature  of  a  real  blade 
decreases  sharply  toward  the  leading  edge,  and  in  a  rigorous  analysis, 
this  should  be  considered.  In  the  present  approach,  however,  the  blade 
surface  has  been  assumed  to  have  a  simply  curved  surface  from  a  leading 
edge  at  x  -  0. 

In  these  analyses,  the  application  of  pressure  gradients  is  independent  of 
surface  curvature,  whereas  aerodynaraically  this  is  not  the  case.  Properly 
sophisticated  usage  of  these  methods  would  entail  supplying  detailed 
pressure  distributions  in  conjunction  with  the  corresponding  surface  raulii 
of  curvature.  The  present  simple  approach  is  employed  to  separate  the 
effects  of  these  various  conditions. 

U.U  Effect  of  Pressure  Gradients 


Pressure  gradients  of  the  form  described  in  Section  2.3.1  were  applied  to 
the  large  and  small  blades  to  determine  the  behavior  of  a  boundary  layer 
on  a  rotor  blade  under  typical  operating  conditions.  The  effects  on  the 
displacement  thickness  are  depicted  in  Figures  (7)  and  (13),  the  corre¬ 
sponding  skew  angle  plots  are  in  Figures  (19)  and  (25),  and  the  velocity 
profiles  are  in  Figures  ( 31 )  and  (37)-  Relevant  data  are  summarized  in 
Tables  III  and  IX. 

Over  the  first  quarter  of  the  chord,  the  free  stream  velocity  is  uniformly 
higher  than  it  is  for  the  zero  pressure  gradient  cases.  This  causes  a 
reduction  ir.  the  boundary  layer  thickness  and  in  the  shear  stress,  as  is 
expected.  At  the  quarter  chord,  the  velocity  gradient  region  begins  and 
the  boundary  layer  is  subjected  to  a  decelerating  external  flow.  The 
increasing  pressure  acts  in  the  same  direction  as  the  shear  stress,  further 
retarding  the  flow  and  causing  an  increased  growth  rate  of  the  boundary 
layer.  There  is  an  associated  lowering  of  the  normal  velocity  gradients 
and  thus  an  increase  in  the  shear  stress  near  the  surface.  The  boundary 
layer  continues  to  grow  rapidly,  but  the  momentum  of  the  viscous  flow 
increases  more  gradued.lv.  A  relative  measure  of  these  two  quantities  is 
the  shape  factor,  H  =  6^/0xx*  In  the  integral  method,  the  shape  factor 
controls  the  growtn  of  the  boundary  layer  via  the  empirical  von  Doenhoff- 
Tetervin  relation.  In  this  type  of  method  where  the  exact  velocity  profile 
is  unknown,  li  is  used  to  evaluate  the  location  of  the  separation  point  mak¬ 
ing  a  comparison  with  analogous  experimental  results .  When  a  predetermined 
value  is  exceeded,  the  flow  is  said  to  have  separated.  If  the  velocity 
profiles  are  known,  as  in  the  differential  method,  the  location  of  separa¬ 
tion  is  det.  --mined  by  the  occurrence  of  reverse  flow  near  the  surface,  or 
more  precisely,  when  the  gradient  of  the  chordwise  velocity  normal  to  the 
surface  it  the  wall  passes  through  zero. 
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From  Figures  (7)  and  (13),  the  boundary  layer  growth  as  represented  by  6X 
is  seen  to  be  very  rapid  after  the  quarter  chord.  For  the  integral  method, 
separation  occurs  in  all  cases  before  the  trailing  edge.  The  thickness  6 
grows  in  a  similar  manner  while  the  momentum  thickness  9xx  increases  less 
rapidly.  A  nonrotating,  two-dimensional  boundary  layer  with  a  free  stream 
velocity  which  is  identical  to  the  corresponding  chordwise  velocity  at 
any  spanwise  location  has  been  found  to  exhibit  the  same  characteristics ; 
rotaticr  does  not  appear  to  affect  the  growth  and  separation  of  the  viscous 
flow  in  the  integral  method. 

In  Figure  ( 1*5 ) ,  the  growth  of  the  shape  factor  is  presented  as  a  function 
of  the  chordwise  position  for  both  pressure  gradient  cases  on  the  large 
blade.  As  the  shape  factor  nears  a  value  of  2.0  for  the  integral  method, 
its  gradient  increases  very  rapidly.  The  value  of  2.0  was  assumed  to  be 
the  separation  point  for  this  method.  This  is  within  the  range  of  H  =  1.8 
to  H  =  2.6  suggested  by  von  Doenhoff  and  Tetervin  for  the  occurrence  of 
separation.  The  slope  of  H  in  this  range  is  such  that  the  location  of  the 
separation  point  is  not  significantly  affected  by  the  value  chosen  as  the 
criterion  for  separation.  In  the  spanwise  direction,  H  was  found  to  in¬ 
crease  toward  the  tip  in  the  integral  method  when  the  adverse  pressure 
gradients  were  applied.  This  is  in  part  related  to  the  increase  in 
caused  by  compressibility,  but  more  importantly,  it  is  a  function  of  the 
rate  of  deceleration  of  the  flow.  Although  the  pressure  coefficients  are 
constant  in  the  spanwice  direction,  the  increasing  rotational  velocity 
yields  stronger  adverse  pressure  gradients,  and  therefore  velocity 
gradients,  toward  the  tip.  This  causes  the  increased  growth  rate  and 
earlier  separation  compared  to  inboard  regions. 

The  skew  angles  are  noticeably  smaller  in  the  pressure  gradient  cases,  even 
over  the  constant  pressure  region  ahead  of  the  quarter-chord.  It  should  be 
remembered  that,  for  the  integral  method,  e  was  found  to  be  an  inverse 
function  of  the  spanwise  position  y  and  relatively  independent  of  free 
stream  velocity  or  angular  velocity.  The  reduction  in  e  is  due  to  the 
increase  in  U»  relative  to  fiy.  In  Equation  (l7y),  this  has  the  effects  of 
increasing  the  Coriolis  term  while  the  centrifugal  term  remains  unchanged. 
This  reduces  3v/3x,  which  in  turn  results  in  a  lower  skew  angle. 


The  shape  factor  variation,  however,  does  not  follow  that  of  the  integral 
method;  a  comparison  of  the  two  methods  is  given  in  Figure  (^5).  The  zero 
pressure  gradient  value  calculated  by  the  differential  method  is  slightly 
higher  them  the  constant  1.286  assumed  in  the  integral  method.  The  calcu¬ 
lated  value,  however,  remains  constant  over  the  chord.  When  pressure 
gradient  is  present,  the  differential  method  does  not  show  the  rapid  diver 
gence  predicted  by  the  integral  method.  If  shape  factor  is  to  be  used  as 
an  indicator  of  the  proximity  to  separation,  then  the  variation  of  H  pre¬ 
dicted  by  the  differential  method  is  consistent  with  the  variation  in  skin 
friction  calculated,  since  separation  is  not  indicated  in  any  of  the  cases 
Figure  (U6)  shows  the  comparable  surface  shear  stress  variations  for  both 
the  integral  and  differential  methods.  In  view  of  the  comparatively  good 
agreement  between  the  results  from  the  two  methods ,  it  may  be  concluded 
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that  the  shape  factor  form  used  in  the  integral  method  is  inappropriate  in 
the  three-dimensional  flows  examined. 

No  experimental  data  is  available  to  correlate  the  separation  position  for 
the  rotating  case,  but  unpublished  two-dimensional  results  indicate  that 
some  separation  should  occur  for  the  3Cp/3(x/c)  =  2.0  cases  close  to  the 
trailing  edge.  The  fact  that  it  dpes  not  occur  in  the  calculated  cases, 
but  is  obviooly  close,  (Figure  (46)),  is  an  indication  of  the  pseudo 
favorable  pressure  gradient  effect  introduced  by  the  radial  gradients  of 
chordwise  velocity  and  the  chordwise  gradients  of  spanwise  velocity. 

4.5  Combined  Effects  of  Pressure  Gradient, 

Crossflow,  and  Curvature 

The  final  effort  in  the  present  study  is  a  complete  rotor  blade  simulation 
for  which  the  effects  of  surface  curvature,  pressure  gradient,  and  circu¬ 
lation  strengths  were  used  on  both  blades  with  the  intermediate  angular 
velocities.  These  conditions  are  felt  to  be  a  reasonable  duplication  of 
a  typical  rotor.  The  resultant  displacement  thicknesses  are  shown  in 
Figures  (10J  and  (l6),  the  skew  angles  sore  plotted  in  Figures  (22)  and 
(20),  and  the  velocity  profiles  are  shown  in  Figures  (34)  and  (40).  The 
summary  data  is  tabulated  in  Tables  VI  and  XII. 

The  influence  of  surface  curvature  is  negligibly  small  for  both  the  vortex 
crossflow  and  the  pressure  gradient  cases.  Its  effect  is  independent  of 
either  of  these  two  conditions,  and  thus  the  combination  of  them  is  not 
expected  to  induce  any  unusual  surface  curvature  effects.  Therefore,  only 
the  relationship  of  the  results  to  the  applied  pressure  gradients  and 
vortex  crossflows  will  be  considered. 

On  the  inboard  portion  of  the  blade  where  the  influence  of  the  vortex  is 
small,  the  results  are  substantially  the  same  as  the  previous  pressure 
gradient  cases.  The  outer  portion  experiences  the  full  strength  of  both 
flow  conditions.  The  effect  on  the  spanwise  distribution  of  the  displace¬ 
ment  thickness  and  the  skew  angle  is  shown  in  Figure  (47).  These  plots 
are  similar  to  Figure  (43)  in  which  only  the  crossflow  was  present,  but 
the  change  to  am  optimum  scale  size  should  be  noted,  as  well  as  the 
difference  in  chordwise  position.  Separation  forced  the  use  of  an  earlier 
position  in  the  present  case.  Nevertheless,  a  comparison  of  the  trends 
reveals  distinctive  features  for  each  situation. 

The  displacement  thickness  increases  in  the  spanwise  direction  with  the 
pressure  gradient  condition;  this  compressibility  effect  is  discussed  in 
Sections  4.1  and  4.3.  The  presence  of  the  interfering  vortex  causes  eui 
alternate  thinning  and  thickening  of  the  boundary  layer  as  it  did  without 
the  pressure  gradient.  The  magnitude  of  the  effect  is  about  the  same  in 
either  case,  but  the  proportion  of  the  thickness  deviation  is  less  in  the 
combined  vortex,  pressure  gradient  case.  This  is  caused  by  the  higher 
chordwise  velocity  which  reduces  the  significance  of  the  imposed  spanwise 
velocity.  The  boundary  layer  thickness  and  the  momentum  thickness  both 
decrease  with  span  because  they  are  relatively  unaffected  by  compressibil¬ 
ity;  in  this  respect  they  are  more  similar  to  the  zero  pressure  gradient 
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cases,  although  the  general  thickness  level  is  greater. 

The  skev  angle  for  the  integral  method  exhibits  a  remarkable  change  between 
the  previous  vortex  cases  and  the  present  combined  cases  .  The  relatively 
small  inflection  noted  then  becomes  masked  by  a  large  skew  angle  with  a 
peak  value  near  the  centerline  of  the  vortex  and  a  sign  which  is  the  same 
as  that  of  the  imposed  external  spanwise  flow.  This  behavior  arises  from 
a  term  containing  the  product  Uoo(3Voo/3y)  in  the  complete  expansion  of  the 
spanwise  integral  equation.  Its  occurrence  is  related  to  the  form  of  the 
power  law  velocity  profiles;  its  practiced  significance  is  dubious. 

Skew  angle  variation  predicted  for  these  flows  by  the  differential  method 
is  largely  the  same  a3  that  for  the  zero  pressure  gradient,  vortex  cross- 
flow  cases  discussed  earlier;  the  form  of  the  curves  is  identical,  and  the 
scale  only  varied  in  chordwise  external  velocity  and  surface  shear  stress. 


The  chordwise  shear  stress  retains  its  sensitivity  to  the  crossflow  noted 
above,  and  vortex  induced  reduction  brings  the  flow  even  closer  to  separa¬ 
tion  than  the  chordviae  pressure  gradient  indicates.  This  is  highlighted 
in  Figure  (44),  vhere  the  possibility  exists  that  an  interfering  vortex 
could  cause  local  separation  if  the  chordwise  pressure  were  any  larger  than 
the  value  of  3Cp/3(x/c)  *  2.0  used  here. 


The  Velocity  Profiles 


In  addition  to  providing  overall  boundary  layer  development,  the  differ¬ 
ential  method  printed  out  at  selected  points  the  velocity  profiles  as  a 
function  of  both  z  and  z/6.  In  addition  to  the  chordwise  and  spanwise 
velocity  ratios,  c  and  s,  presented  in  Figures  (29)  through  (4o),  stream- 
wise  and  normal  velocity  profiles  were  available.  There  is  nothing  remark¬ 
able  about  the  chordwise  and  spanwise  velocity  profiles  presented.  The 
chordwise  profile  follows  a  typical  two-dimensioned  form  closely  matching 
the  typical  1/7  power  law  form  commonly  assumed.  The  spanwise  profile,  as 
far  as  con  be  determined  from  the  data  available,  follows  the  external 
spanvise  velocity,  with  a  form  closely  resembling  the  chordwise  profile. 

The  parallels  between  the  rotating  and  two-dimensional  velocity  profiles 
are  reinforced  by  Figure  (48)  where  a  typical  chordwise  profile  is  compared 
with  the  standard  law  of  the  wall  form.  Agreement  is  remarkably  close, 
the  effect  of  the  rotation  being  essentially  negligible.  When  a  pressure 
gradient  case  is  compeared,  the  characteristic  wake  contribution  appears  in 
the  outer  levels  of  the  boundary  layer.  The  addition  of  a  large  external 
crossflow,  vortex  induced,  causes  only  minor  changes  in  the  wake  portion  of 
the  flow.  For  Figure  (48),  the  chordwise  shear  velocity  was  defined  as 
UT  *  (tx/Pod)**.  The  effect  of  the  pressure  gradient  on  the  velocity  profile 
in  the  typical  z/6  vs  u/U»  plot  is  clear  in  Figure  (49)  with  the  increased 
fullness  of  a  profile  approaching  separation. 


Velocity  profiles  were  printed  out  at  x/chord  =  0.3  and  0.8;  consequently, 
no  profiles  of  flows  at  or  immediately  prior  to  separation  are  presented. 
Figure  (50)  gives  the  typical  variation  of  the  chordwise,  spanwise,  and 
normal  velocity  ratios  for  a  case  with  zero  pressure  gradient  and  no  vortex 
crossflow.  The  normal  velocity  component  does  not  follow  directly  the 


external  crossflow.  At  x/ chord  =  0.3,  the  spanwise  ratio,  s,  has  already 
crossed  over,  following  the  changes  in  V»,  and  is  negative.  The  normal 
velocity,  n,  is  still  positive.  By  the  time  x/chord  =  0.8,  n  has  also 
become  negative.  This  tendency  for  the  normal  component  sign  to  lag  the 
external  flow  is  more  marked  for  the  inboard  sections  and  is  particularly 
noticeable  on  the  small  blade. 

The  general  for?>  of  the  profiles,  if  not  the  sign  of  the  crossflow  profile, 
agrees  with  those  assumed  for  the  integral  method  and  presented  in  Figure 

(2). 

4 .  7  The  Shear  Stress  Profiles 

As  with  the  velocity  profiles,  there  was  very  little  to  distinguish  the 
chordwise  shear  stress  profiles  from  similar  profiles  in  an  equivalent 
nonrotating  flow.  Figure  (51)  shows  typical  chordwise  shear  stress  pro¬ 
files  for  two  cases  with  and  without  pressure  gradient.  The  zero  pressure 
profile  has  the  characteristic  triangular  shape,  with  the  shear  stress 
falling  away  from  the  edge  of  the  laminar  sublayer.  When  the  flow  is  sub¬ 
jected  to  an  unfavorable  pressure  gradient,  the  shear  stress  profile  takes 
on  the  characteristic  form  of  Figure  (50).  The  shear  stress  increases 
away  from  the  wall  to  a  maximum  of  about  twice  the  wall  value  at  z/6  =  0.25. 
Outside  this,  the  shear  stress  falls  off  to  zero  in  the  external  stream. 
Since  the  shear  stresses  in  the  turbulent  part  of  the  layer  are  determined 
by  resolving  the  total  shear  stress  based  on  the  resultant  flow,  the  span- 
wise  and  chordwise  shear  stresses  are  related  through  the  ratios  3s/3z  and 
3c/3z.  The  spanwise  stress  normal  variation  is  consequently  similar  to 
that  of  the  chordwise  stress;  the  two  values  away  from  the  wall  are  of 
relative  magnitude  Cfy/cp*. 

It  is  the  gradient  of  shear  stress  normal  to  the  wall,  rather  than  the 
shear  stresses  themselves,  that  determines  the  boundary  layer  growth.  The 
variation  of  the  shear  stress  gradient  normal  to  the  wall  for  the  two  cases 
discussed  above  is  presented  in  Figure  (52).  The  portions  of  these  curves 
outside  the  sublayer  have  been  obtained  by  differentiating  the  shear  stress 
profiles.  The  balance  of  the  curve  is  found  from  the  interpolation  based 
on  u+2  discussed  in  Section  3.1.1.  The  portion  of  the  curve  close  to  the 
wall  shown  in  an  expanded  scale  illustrates  the  importance  of  the  form  of 
the  interpolation,  particularly  in  the  pressure  gradient  case  where  large 
positive  values  of  3t/3 l  are  present  close  to  the  wall. 

4.8  General  Discussion 

The  results  predicted  for  the  characteristic  thickness  parameters  and  the 
chordwise  surface  shear  stress  coefficient  were  in  remarkably  good  agree¬ 
ment  for  all  the  cases  where  there  was  zero  chordwise  pressure  gradient 
There  was,  however,  a  divergence  in  the  predicted  values  of  skew  angle. 

Both  methods  predicted  an  initially  positive  value  ahead  of  the  quarter 
chord  line.  The  differential  method  showed  a  tendency  toward  negative 
values  behind  this  point,  whereas  the  integral  method  continued  to  show 
growing  positive  values.  A  typical  example  is  in  Case  2  with  0  *  15 
rad/sec,  3Cp/3(x/C)  =  0  and  T  =  0.  At  x/chord  =  1.0  and  y/R  ■  0.3,  the 
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integral  and  differential  methods  show  e  =  0.590  and  -  0.255  degrees  re¬ 
spectively.  The  apparent  difference  is  smaller  when  it  is  realized  that 
these  represent  surface  flow  angularity  of  -6.536  and  -7.381  degrees  rela¬ 
tive  to  the  free  stream  velocity  in  a  direction  of  -7.126  degrees.  The 
correlation  between  the  surface  shear  stresses  in  the  spanwise  direction 
is  fair  with  values  of  -0.000229  and  -0.000183  from  the  differential  and 
integral  methods. 

The  small  difference  between  the  skew  angles  predicted  by  the  two  methods 
is  basic  and  results  from  the  assumption  in  the  integral  method  that  the 
crossflow  velocity  profile  is  coupled  to  the  streamvise  flow  through  a 
function  containing  s.  The  questions  arising  from  this  assumption  are 
fully  discussed  in  Reference  (9).  They  must  be  kept  in  mind  when  evalu¬ 
ating  the  skew  angles  predicted  by  the  integral  method. 

The  difference  in  predicted  skew  angle  is  obvious  in  Figures  (U3)  and  (U7) 
where  the  effect  of  vortex  crossflow  is  examined  with  and  without  the  in¬ 
fluence  of  chordwise  pressure  gradient.  The  inability  of  the  integral 
method  to  accurately  predict  skew  angle  is  highlighted  here  with  an 
apparently  significant  influence  from  imposed  pressure  gradient.  It  is 
felt,  on  the  basis  of  the  above  discussion  and  the  similar  form  of  the  skew 
angles  predicted  by  the  differential  method  with  and  without  pressure 
gradient,  that  the  differential  method  givus  the  more  appropriate  descrip¬ 
tion  of  what  is  happening. 

The  two  methods  also  differ  on  the  position  at  which  the  chordwise  surface 
shear  stress  becomes  zero.  The  integral  method  invokes  the  assumption  that 
when  the  shape  factor  H  equals  two,  then  the  flow  has  separated.  This  is 
inappropriate  in  three-dimensional,  rotating  flow,  as  seen  in  Figures  (1+5) 
and  (U6)  where  the  chordwise  variations  of  H  and  cfX  are  plotted  for  a 
range  of  pressure  gradients  for  the  two  methods.  Even  after  a  value  of  H 
greater  than  2.0  has  been  reached  in  the  integral  method  and  the  flow  has 
theoretically  separated,  the  correlation  between  the  skin  friction  values 
predicted  by  the  two  methods  is  still  fair.  On  the  oasis  of  this,  the 
shape  factor  skin  friction  relationships  used  are  suspect  in  the  flows 
presented  here  and  should  be  reexamined  before  extending  the  study  to  more 
involved  pressure  distributions .  They  appear  to  be  adequate  in  the  zero 
chordwise  pressure  gradient  cases. 

No  reliable  experimental  data  is  available  to  correlate  the  separation 
points  calculated,  although  unpublished  two-dimensional  data  indicates 
that  separation  might  be  expected  close  to  the  trailing  edge  for  the 
3Cp/3(x/C)  =  2.0  cases  on  the  large  blade.  That  it  is  not  predicted  using 
the  differential  method  could  be  an  alleviation  of  the  unfavorable  pressure 
gradient  by  the  effects  of  the  rotation,  or  a  failure  of  the  method.  A 
final  decision  on  this  point  cannot  be  made  until  more  realistic  pressure 
distributions  are  calculated  and  relevant  experimental  data  becomes 
available. 


CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 


The  present  analyses  indicate  that  on  rotating  flat  or  slightly  curved 
surfaces  of  typical  helicopter  rotor  size,  in  the  absence  of  chordwise 
pressure  gradients,  the  boundary  layer  behaves  in  an  essentially  two- 
dimensional  manner  when  followed  along  the  streamlines  of  the  external 
flow. 

The  influence  of  large  amounts  of  vortex  induced  crossflow  on  the  zero 
pressure  gradient  rotating  boundary  layer  is  small. 

Centrifugal  force  has  very  little  influence  on  the  flow  at  the  angular 
velocities  examined  and  only  becomes  moderately  significant  when  compress¬ 
ibility  is  a  controlling  factor  at  high  subsonic  Mach  numbers. 

Coriolis  force  is  the  most  significant  force  present.  Together  with  the 
radial  variation  of  chordwise  velocity,  it  dominates  the  development  of 
the  spanwise  flow. 

The  influence  of  chordwise  pressure  gradients  is  reduced  by  the  presence 
of  a  negative  spanwi.se  external  flow.  This  effect  is  most  noticeable  at 
small  radii  and  close  to  the  trailing  edge  where  the  spanwise  flow  is 
large  relative  to  the  chordwise  velocity. 

Both  the  differential  and  integral  methods  agree  remarkably  well  on  the 
prediction  of  the  boundary  layer  characteristic  thicknesses  for  the  range 
of  conditions  examined.  Correlation  between  the  methods  for  the  skin 
friction  values,  in  all  but  the  pressure  gradient  cases,  is  also  good.  It 
may  be  concluded  that  the  empirical  skin  friction  lavs  used  in  the  integral 
technique  are  applicable  in  three-dimensional  flows,  at  least  for  zero 
pressure  gradient. 

It  would  be  inappropriate  to  conclude  that  since  the  two  widely  differing 
methods  give  results  which  are  generally  in  good  agreement  that  they  are 
capable  of  giving  a  complete  picture  of  turbulent  boundary  layers  on 
helicopter  rotors.  It  is  however  encouraging  that  such  good  agreement  has 
been  reached,  over  a  wide  range  of  conditions,  using  two  comparatively 
basic  analyses  with  the  minimum  of  restrictive  assumptions. 

Despite  the  question  marks  which  must  remain  over  the  analyses  until  con¬ 
firmatory  experimental  data  becomes  available ,  it  should  now  be  possible  to 
more  realistically  allow  for  the  effects  of  rotation  on  airfoil  skin 
friction  and  maximum  lift  capability. 

RECOMMENDATIONS 


Before  extending  the  current  analyses  to  unsteady  flows,  or  even  more 
exotic  steady  flows,  some  experimental  data  on  boundary  layer  development 
on  rotating  blades  must  be  found.  Data  needed  include  velocity  profiles, 
surface  shear  stresses  and  thickness  parameters,  as  well  as  some  measurement 
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of  the  fluctuating  velocities  and  the  determination  of  turbulent  shear 
stress  variation  normal  to  the  wall  in  rotating  flows .  These  data  are 
essential  for  meaningful  analysis  of  the  present  analytical  work  and  any 
future  developments  of  the  theory. 

In  addition,  the  use  of  the  present  methods  can  be  expanded  by  the  follow¬ 
ing  steps.  Extend  the  analysis  from  its  present  form  now  handling  simply 
curved  surfaces  to  one  capable  of  handling  the  sharp  leading  edge  curves 
of  airfoils.  Develop  the  capability  of  calculating  complete  viscous  flows 
from  stagnation  point,  through  laminar  and  transition  phases,  to  the  fully 
developed  turbulent  flow  in  a  rotating  field.  Modify  the  present  analyses 
to  allow  the  calculation  of  the  behavior  of  the  boundary  layer  in  unsteady 
conditions  and  particularly  in  the  forward  flight  regime. 
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X  -  POSITIVE  MEASURED 
FROM  QUARTER  CHORD 
PARALLEL  TO  SURFACE 


Figure  1.  Rotating  Orthogonal  Curvilinear 
Coordinate  System. 
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PRESSURE  COEFFICI 


-1.0  H 


APPROXIMATED  PRESSURE 

DISTRIBUTION  ,  4%  =  ' 

’  dx/C 

MEASURED  PRESSURE  DISTRIBUTION  , 
ANGLE  OF  ATTACK  =  0°  , 

MACH  NUMBER  =  0.65 


Figure  3.  Comparison  of  Approximated  Pressure  Distribution 
end  Measured  H-3l*  Airfoil  Data. 
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CURVED  APPROXIMATION 
NACAOOI2  SURFACE 
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DISPLACEMENT  THICKNESS, 


.008 


INTEGRAL  METHOD 


t  >  1-1  A  I  I  1  -  A  i 

0  .2  .4  .6  .8  1.0  (0.95) 

CHORDWISE  POSITION,  x/C 


(b)  ft  =  15  rad/sec 
Figure  5 •  Continued. 
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DISPLACEMENT 


*  »  *  *  *  *  »  *  *  *  i 

0  .2  .4  .6  8  1.0  (0.95 

CHORDWISE  POSITION,  x/C 


(a)  r/h  =  +100  ft/sec 

Figure  6.  Displacement  Thickness  Versus  Chordwise  Position  for 
Various  Vortex  Strengths,  Blade  Radius  =  Lo  ft, 

=  15  rad/sec,  dCp/d(x/C)  =  0,  Flat  Surface. 
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DISPLACEMENT  THICKNESS,  8 
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INTEGRAL  N  .THOD 


^  *  *  1  *  *  *  *  A  ^  » 

0  .2  .4  .6  .8  1.0  |jQ.95) 

CHORDWISE  POSITION,  x/C 


(c)  r/h  =  +200  ft/sec 
Figure  6.  Continued. 
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INTEGRAL  METHOD 


(d)  r/h  =  -200  ft/sec 
Figure  6.  Concluded. 
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CHORDWISE  POSITION,  x/C 


(a)  dCp/d(x/C)  =  1 

Figure  7.  Displacement  Thickness  Versus  Chordwise  Position  for 
Various  Pressure  Gradients,  Blade  Radius  =  ^0  ft, 

0  =  15  rad/sec,  r/h  =  0  ft/sec.  Flat  Surface. 


95 


DISPLACEMENT  THICKNESS,  8 
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0  .2  .4  .8  .8  1.0  (13.95) 

CHORDWISE  POSITION,  x/C 

(a)  T/h  =  +100  ft/sec 

Figure  8.  Displacement  Thickness  Versus  Chordvise  Position  for 
Various  Vortex  Strengths,  blade  Radius  =  Uo  ft, 
ft  =  15  rad/sec,  dCp/d(x/C)  =  0,  Surface  Radius  of 
Curvature  =  8.25  ft. 
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CHORDWISE  POSITION,  x/C 


(b)  r/h  =  -100  ft/sec 
Figure  8.  Continued. 


98 


(c)  T/h  =  +200  ft/sec 
Figure  8.  Continued. 


Figure  8.  Concluded. 


DISPLACEMENT  THICKNESS,  3 


-  .004 

*x 


INTEGRAL  METHOD 


SEPARATION 
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DIFFERENTIAL  METHOD 


CHORDWISE  POSITION,  x/C 


(r/R  =  0.30 
(0.8 
to  fc).9 
"to  (p.9 


(a)  dCp/d(x/C)  =  1 

Figure  9*  Displacement  Thickness  Versus  Chordvise  Position  for 
Various  Pressure  Gradients,  Blade  Radius  =  1*0  ft, 
n  =  15  rad/sec,  T/h  =  0  ft/sec.  Surface  Radius  of 
Curvature  =  8.25  ft. 


DISPLACEMENT  THICKNESS,  8 


.016 


(b)  dCp/d(x/C)  =  2 
Figure  9*  Concluded. 
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0  .2  4  .8  .8  10  J3.95 

CHORDWISE  POSITION,  x/C 


(a)  r/h  =  +200  ft/sec 

Figure  10.  Displacement  Thickness  Versus  Chordvise  Position  for 
Various  Vortex  Strengths,  Blade  Radius  =  hO  ft, 
ft  =  15  rad/sec,  dCp/d(x/C)  =  2,  Surface  Radius  of 
Curvature  =  8.25  ft. 
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6  ‘  t  ‘""4  .6  ‘  .8  1.0  (0.95) 

CHORDWISE  POSITION,  x/C 


(a)  n  =  LO  rad/sec 

Figure  11.  Displacement  Thickness  Versus  Chordwise  Position  for 
Various  Angular  Velocities,  Blade  Radius  =  10  ft, 
f/h  =  0  ft/sec,  dCp/d(x/C)  =  0,  Flat  Surface. 
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DISPLACEMENT  THICKNESS,  8 


DISPLACEMENT  THICKNESS,  8 


CHORDWISE  POSITION,  x/C 


(c)  n  =  80  rad/sec 


Figure  11.  Concluded. 
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CHORDWISE  POSITION,  x/C 

(a)  r/h  =  ►25  ft/sec 

Figure  12.  Displacement  Tnickness  Versus  Chordwise  Position  for 
Various  Vortex  Strengths,  Blade  Radius  =  10  ft, 

Q  =  60  rad/sec,  dCp/d(x/C)  =  0,  Flat  Surface. 
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t  I  *  l  .1  ,  ■  -1  i  I  l  I  l  . 

0  2  .4  .6  .8  1.0  (0.95) 

CHORDWISE  POSITION,  x/C 

(c)  r/h  =  +50  ft/sec 
Figure  12.  Continued. 
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(d)  r/h  =  -50  ft/sec 


Figure  12.  Concluded. 
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CHORDWISE  POSITION,  x/C 

(a)  dCp./d(x/C)  =  1 

Figure  13.  Displacement  Thickness  Versus  Chordwise  Position  for 
Various  Pressure  Gradients,  Blade  Radius  =  10  ft, 
ft  =  60  rad/sec,  f/h  =  0  ft/sec.  Flat  Surface. 
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0  .2  .4  .6  .8  1.0  (0.95) 

CHORDWISE  POSITION,  x/C 

(b)  dCp/d(x/C)  =  2 
Figure  13.  Concluded. 
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^  c )  r  h  =  +50  ft/sec 


Figure  1^-.  Cortinued. 
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Figure  lU.  Concluded. 
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CHORDWISE  POSITiuN,  x/C 

(a)  dCp/d(x/C )  =  1 

Figure  15-  Displacement  Thickness  Versus  Chordvise  Position  for 
Various  Pressure  Gradients,  Blade  Radius  =  10  ft, 

=  60  rad/sec,  r/h  =  0  ft/sec.  Surface  Radius  of 
Curvature  =  U.125  ft. 
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(b)  dCp/d( x/C)  =  2 


Figure  15-  Concluded. 
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DISPLACEMENT  TH 


(b)  r/h  =  -50  ft/sec 
Figure  16.  Concluded. 
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INTEGRAL  METHOD 
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(a)  il  *  10  rad/sec 

Figure  17-  Skew  Angle  Versus  Chordwise  Position  for  Various 
Angular  Velocities,  Blade  Radius  *  Uo  ft, 

T/h  *  0  ft/sec,  dCp/d(x/C)  =  0,  Flat  Surface. 
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(b)  Q  »  15  r*d/sec 
Figure  17-  Continued. 
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(c)  J)  =  20  rad/sec 
Figure  17.  Concluded. 
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SKEW  ANGLE,  €,  DEG 


CHORDWISE  POSITION,  x/C 

(a)  r/h  =  +100  ft/sec 

Figure  18.  Skew  Angle  Versus  Chordwise  Position  for  Various 
Vortex  Strengths,  Blade  Radius  =  ho  ft, 
n  =  15  rad/sec,  dCp/d(x/C)  =  0,  Flat  Surface. 
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0  .2  .4  .6  .8  1.0  (0.95} 

CHORDWISE  POSITION,  x/C 


(b)  r/h  =  -100  ft/sec 
Figure  18.  Continued. 
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Figure  18.  Continued. 
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Figure  18.  Concluded. 
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Figure  19.  Skew  Angle  Vereu*  Chordwise  Position  for  Various 
Pressure  Gradients,  Blade  Radius  ■  UO  ft, 

Q  ■  15  rad/sec,  f/h  ■  0  ft/sec.  Flat  Surface. 
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CHORDWISE  POSITION,  x/C 


(b)  dCp/d(x/C)  =  2 
Figure  19.  Concluded. 
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SKEW  ANGLE,  €,  DEG 


SKEW  ANGLE,  €,  DEG 


(b)  r/h  =  -100  ft/sec 
Figure  20.  Continued. 
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Figure  21.  Skew  Angle  Versus  Chordwise  Position  for  Various 
Pressure  Gradients,  Blade  Radius  =  10  ft , 
ft  =  15  rad/sec,  T/h  =  0  ft/sec.  Surface  Radius 
of  Curvature  =  8.25  ft. 
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Figure  21.  Concluded. 
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Figure  22.  Skew  Angle  Versus  Chordwise  Position  for  Various 
Vortex  Strengths,  Blade  Radius  =  Uo  ft, 
ft  =  15  rad/sec,  dCp/d(x/C)  =  2,  Surface  Radius 
of  Curvature  =  8.25  ft. 
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Skew  Angle  Versus  Chordwise  Position  for  Various 
Angular  Velocities,  Blade  Radius  =  10  ft, 
f/h  =  0  ft/sec,  dCp/d(x/C)  =  0,  Flat  Surface. 
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(b)  dCp/d(x/c)  =  2 
Figure  25.  Concluded. 
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Vortex  Strengths,  Blade  Radius  =  10  ft, 
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Figure  Continued. 
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Figure  27.  Skew  Angle  Versus  Chordwise  Position  for  Various 
Pressure  Gradients,  Blade  Radius  =  10  ft, 

0  =  60  rad/sec,  T/h  =  0  ft/sec.  Surface  Radius 
of  Curvature  =  U.125  ft. 
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(a)  r/h  =  +50  ft/sec 

Figure  28.  Skew  Angle  Versus  Chordwise  Position  for  Various 
Vortex  Strengths,  Blade  Radius  =  10  ft, 

Q  =  60  rad/sec,  dCp/d(x/C)  =  2,  Surface  Radius 
of  Curvature  =  4.125  ft. 
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Figure  29. 
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Figure  29 •  Continued. 
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Figure  29.  Continued. 
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Figure  29.  Concluded. 
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Figure  30.  Continued. 
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Figure  31.  Continued. 
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Fig 'ore  3<- .  Chordvise  and  cpar.vise  Velocity  Profiles  for  Various  Vortex 
Strengths,  Blade  Fadius  =  1*0  ft,  ft  =  15  rad/sec, 
dCpd(x'C)  =  0,  Surface  Radius  of  Curvature  =  8.25  ft. 
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Figure  32.  Continued. 
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Figure  32.  2ontir.uei. 
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Figure  33.  Chordwise  and  Spanwise  Velocity  Profiles  for  Various  Pressure 
Gradients,  Blade  Radius  =  1+0  ft,  fl  =  15  rad/sec, 

T/h  =  0  ft/sec.  Surface  Radius  of  Curvature  =  8.25  ft. 
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Figure  33-  Continued, 


x/C=0.30 

r/R=0.90 
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unorawise  ana  bpanwise  velocity  Profiles  for  Various  Vortex 
Strengths,  Blade  Radius  =  Uo  ft,  ft  =  15  rad/sec, 
dCp/d( x/C )  =  2,  Surface  Radius  of  Curvature  =  8.25  ft. 


Figure 


x/C=0.30  I  x/C=0.80 


8/Z  '33VJanS  3A08V  1H3I3H 


190 


Figure  35-  Chordwise  and  Spanwise  Velocity  Profiles  for  Various  Angular 
Velocities,  Blade  Radius  =  10  ft,  T/h  =  0  ft/sec, 
dCp/d(x/C)  =  0,  Flat  Surface. 
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Figure  35 •  Continued. 
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Figure  35.  Continued. 
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Figure  35.  Continued. 
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Figure  36.  Continued. 
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Figure  36.  Continued. 
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Figure  36.  Continued. 
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Figure  36.  Continued. 
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Figure  3 6.  Continued. 
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Figure  36.  Concluded. 
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Figure  37.  Chorth/lse  end  Span vise  Velocity  Profile*  fbr  Various  Pressure 
Gradients,  Blade  Radius  ■  10  ft,  Q  ■  60  rad/sec, 
f/h  ■  0  ft/sec.  Flat  Surface . 
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Figure  37.  Continued 


Figure  38.  Chordvise  and  Spanvise  Velocity  Profiles  for  Various  Vortex 
Strengths,  Blade  Radius  =  10  ft,  n  =  60  rad/sec, 
dCp/d(x/C)  =  0,  Surface  Radius  of  Curvature  =  U.125  ft. 
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(a)  Concluded. 
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Figure  38.  Continued. 
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Figure  38.  Continued. 
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Figure  39-  Chordvise  and  Spanwise  Velocity  Profiles  for  Various  Pressure 
Gradients,  Blade  Radius  =  10  ft,  fl  =  60  rad/sec, 

T/h  =  0  ft/sec.  Surface  Radius  of  Curvature  =  U.125  ft. 
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Figure  40.  Chordwise  and  Spanwise  Velocity  Profiles  for  Various  Vortex 
Strengths,  Blade  Radius  =  10  ft,  fi  =  60  rad/sec, 
dCp/d(x/C)  =  2,  Surface  Radius  of  Curvature  =  4.125  ft. 
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Figure  Ul.  Typical  Effect  of  Angular  Velocity  on  Skew 
Angle  Predicted  by  the  Two  Methods. 
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Figure  h2.  Effect  of  Angular  Velocity  on  Chordvise  Skin  Friction. 
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Figure  Typical  Effect  of  Pressure  Gradient  on  Chordvise 

Variation  of  Shape  Factor,  H,  9Cp/9(x/C)  =  0,  1,  2, 
y/R  =  0.9,  0  =  15  rad/sec,  f  =  0. 
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Figure  hj.  Typical  Effect  of  Vortex  Crossflow  on  Spanvise 

Variation  of  Displacement  Thickness  and  Skew  Angle 
^Cp/9(x/C)  =  2,  x/ Chord  =  0.7»  fi  =  15  rad/sec. 


SS3NX0IH1  1  N3W3D\ndSIQ 


DISTANCE  NORMAL  TO  SURFACE  , 


1 

[ 

Figure  *49.  Typical  Effect  of  Pressure  Gradient  on  Chordvise 
Velocity  Profiles,  3Cp/3(x/C)  =  0,  2, 
x/Chord  =  0.8,  y/R  =  0.9,  fl  =  15  rai/sec,  r  =  0. 
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Figure  51*  TypicaJ  Variation  of  Chordwise  Shear  Stress 
Normal  to  the  Surface,  x/Chord  =  0.8, 
y/P  ■  0.6,  11  ■  60  rad/sec,  f  =  0. 
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Figure  52.  Typical  Variation  of  the  Normal  Gradients  of 
Chordwise  Shear  Stress  Normal  to  the  Surface. 
x/Chord  =  0.8,  y/R  =  0.6,  fi  =  60  rad/sec,  f  =  0 
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APPENDIX 


Determination  of  Sublayer  Shear  Stress  Gradients 


It  has  been  demonstrated  by  several  authors,  but  most  recently  in 
Reference  (20),  that  in  two-dimensional  flow  the  shear  stress  gradient 
normal  to  the  wall  may  be  represented,  in  the  region  close  to  the  wall,  by 
the  expression 

=  Ao  +  A2U+2 

where  A0  and  A2  are  constants  obtained  from  the  equation  of  motion  eval¬ 
uated  at  the  wall  and  u+  is  the  commonly  employed  velocity  ratio  u/ut.  If 
the  law  of  the  wall  is  assumed  valid  in  three-dimensional  flow,  then 
similar  relations  may  be  derived  in  three-dimensional  flows.  They  will  be 
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Expanding  Equations  (63)  and  substituting  for  the  pressure  gradients  gi /es. 
the  expanded  form 
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Terms  involving  gradients  of  the  chordwise  shear  stress  in  the  spanwise 
direction  and  all  gradients  of  the  spanwise  shear  stress  are  of  small  order 
and  have  been  dropped  from  Equations  (6**). 
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The  equations  may  be  simplified  to 


y ~  =  Aq  +  Apfx(z)  +  A2u+2  (65) 

and 

=  B0  +  B1fy(z)  +  B2u+2 

The  middle  term,  absent  in  the  two-dimensional  form,  is  a  result  of  the 
presence  of  the  Coriolis  and  centrifugal  terms. 

The  fit  to  the  profiles  outside  the  sublayer  was  carried  out  through  the 
u+2  term.  A  curve  was  fitted  through  thr.  five  points  just  outside  the 
sublayer  and  the  coefficients  A2  and  B2  in  Equations  (65)  were  forced  to 
give  agreement.  This  means  essentially  that  the  chordwise  gradient  of  the 
skin  friction  is  being  determined  from  the  normal  gradient  of  shear  stress 
outside  the  sublayer. 
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